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Physicists and the War 


HE importance of physics and the physicist 
in this most mechanical of all wars is be- 
coming daily more obvious. There is evidence 
that the Government are now fully aware of the 
necessities of the situation, and are either taking 
or have in contemplation adequate steps to meet 
it. A certain amount of confusion existed in the 
early days of the war, and in particular the actual 
number of physicists who would be required for 
service either with the forces or in industry was 
grossly underestimated, even by those who were 
in the best position to judge, and a certain num- 
ber of students with respectable but not out- 
standing qualifications were called up and posted 
to units where their training could not be ade- 
quately used. The exigencies of the situation, 
assisted perhaps by representations made by 
Officers of the Institute, and particularly by our 
President, have produced a complete change. All 
men completing a Uni- 


the specialized branches of the Forces are filled. 
I have first-hand knowledge of the efficiency, 
skill and sympathy with which the officials of 
Central Register are carrying out this important 
task. 

The Government are fully alive to the necessity 
of ensuring a continuous flow of trained physi- 
cists into the ‘pool.’ No pressure is at present 
being exercised upon the Technical Committees, 
who have the responsibility for the decision, to 
interrupt a man’s scientific training at an unduly 
early stage, though the stringent scrutiny to 
which each student is now subjected at each 
stage of his training will probably result in the 
partial elimination of the ‘third class Honours’ 
man, for which we may feel no regret. It is under- 
stood that steps are being taken to ensure an 
adequate entry of physics students into the 
Universities from the secondary schools. On the 

other hand, University 





versity degree course 
with physics as one of 
their subjects are inter- 
viewed by a Technical 
Committee, nominated 
by the University or 
University College, 
their capability as- 
sessed, and their prefer- 
ences ascertained. They 
are then transferred to 
a central ‘pool’ from 
which vacancies in in- 
dustry and the Uni- 
versities as well as in 





In these days when more and more 
American physicists are being called 
into defense work, the present situa- 
tion of physicists in Great Britain 
may perhaps be a forecast of things 
to come. With that in mind, the 
accompanying article by Dr. J. A. 
Crowther is reprinted with permis- 
sion from ‘‘Notes and Notices,”” a 
private publication for members of 
the British Institute of Physics. 


departments are to a 
certain extent modify- 
ing and abridging their 
courses to adapt them 
to the urgent needs of 
the times, and the stu- 
dent of to-day and to- 
morrow will be _ less 
widely read than his 
predecessors. In par- 
ticular, he will not 
have had the advant- 
age, which was becom- 
ing fairly general in 
recent years, of com- 
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pleting his studies with a year or so of research 
training. It is very much to be hoped that, after 
the war, the Government will adopt some scheme 
for returning at any rate the best of their recruits 
to the Universities for ‘reconditioning.’ 

It must be emphasized that, while the need 

for more physicists is urgent, the demand at 
present is mainly for junior physicists; for men, 
that is, to carry out highly skilled but neverthe- 
less mainly routine jobs, both in industry and in 
connexion with the fighting forces. The con- 
struction, testing, maintenance and operation of 
new and complicated experimental devices calls 
for knowledge and technical ability which only 
the trained physicist can adequately provide; 
and obviously calls for their employment in large 
numbers. It is certainly not research, in the sense 
in which the word is used in peace time, and it 
is understandable that some of the brilliant 
young physicists called upon to carry it on (some 
of whom would have normally been experiencing 
the thrills and torments of original research) may 
feel that their talents could be employed to 
better purpose. 

A similar feeling undoubtedly exists among the 
younger members of industrial research organiza- 
tions, who, having served a few years apprentice- 
ship in their present jobs would now, except for 
the war, have been reasonably expecting promo- 
tion to more important work. One can sympa- 
thize with them in this temporary halting of 
their career, though it is only fair to point out 
that many young professional men have had 
their careers much more drastically interrupted 
by the demands of military service, and one can 
appreciate still more their sincere desire to devote 
all their talents to the service of the nation. The 
Institute is sparing no effort to ensure that, as 
far as is possible, this desire shall achieve fulfil- 

-ment. It must, however, be remembered that 
this routine work, dull and unimportant as it 
may seem to the individual worker, is in the 








aggregate a vital, possibly the most vital part of 
the war effort; the accuracy and skill with which 
it is carried out will play a decisive part in the 
ultimate event. 

There is a little evidence that some firms, 
harrassed by pressure of work, are not making 
it easy for members of their staffs to seek more 
important posts elsewhere, fearing the difficulties 
of replacement under present conditions. The 
Institute is prepared to advise members, should 
such cases arise; though, where all is so impor- 
tant, questions of ‘priorities’ are not easy to 
decide. The Ministry of Labour possesses wide 
powers in such cases, though it is, quite rightly, 
reluctant to exercise them, and would only act 
on the grounds of national, not personal, ad- 
vantage. In the meantime, we shall continue to 
notify all members whose names are on the 
Appointments Register of vacancies as_ they 
occur. 

Of genuine fundamental war research it is 
unprofitable to pronounce. Very few people can 
be in a position to estimate either its magnitude 
or its quality; and they, for obvious reasons, 
cannot tell. The new, ad hoc, teams have been 
manned by the obvious people from the world of 
pure research, though there is a feeling that a 
larger admixture of men with industrial research 
experience might have expedited the work. There 
may be some justification for this feeling, though 
it has to be remembered that some of our leading 
industrial teams are engaged as such on full-time 
war research. We can, at any rate, welcome the 
fact that the newly appointed Engineering Ad- 
visory Committee has, on its membership, three 
distinguished Fellows of our own Institute, Dr. 
A. P. M. Fleming (recently a member of the 
Board), Dr. C. C. Paterson (ex-President) and 
Sir Henry Tizard. Their presence on the com- 
mittee should ensure that the vital contribution 
which physics still has to offer to the national 
effort will not be disregarded. 
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Résumé of Recent Research 














The Dislocation 
Theory of Plastic 
Deformation 


The first extensive use 
of dislocation theory oc- 
curred in 1934 in the 
papers of Taylor, Oro- 
wan, and Polanyi.! Taylor used the theory to 
discuss work hardening in aluminum, copper, 
and rocksalt. He supposed that, if a crystal is 
subjected to a shearing stress of about 10° 
dynes/cm?*, imperfections of the type shown in 
the figures below are formed somewhere in the 











e 
NEGATIVE DISLOCATION 


Generation of dislocations at crystalline boundaries 
(after Taylor). In cases (a), (b), (c) a positive dislocation 
is generated at the left-hand side of the crystal and moves 
toward the right, whereas in cases (d), (e), (f) a negative 
dislocation moves to the left. The arrows in figure (a) indi- 
cate the component of the external shearing stress which is 
responsible for the plastic deformation. 


crystal and move across the specimen under the 
influence of the external stress. The passage of a 
dislocation entirely across the specimen produces 
a deformation of the material; the upper portion 
of the crystal is displaced relative to the lower 
portion by one atomic spacing. Taylor found 
that like dislocations repel one another whereas 
unlike dislocations attract. Thus as the deforma- 
tion of an initially unstrained crystal proceeds, 
large numbers of dislocations get stuck in the 
interior of the specimen. (The forces acting 
between dislocations do not allow one dislocation 
to pass freely by another dislocation.) A larger 
external stress is then required to move disloca- 
tions through the crystal and the material has 
become work hardened. 
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Seitz and Read, W. F. Brown, Jr., and Mott 
and Nabarro? have examined further conse- 
quences of the theory. Their investigations have 
not revealed any grave failures of the theory. 

In a recent paper J. S. Koehler* has discussed 
additional results which can be derived using 
dislocation theory. It is shown that a dislocation 
near the surface of a specimen experiences an 
“image” force which acts towards the surface. 
Dislocations near the surface will therefore 
diffuse out of a work-hardened specimen and 
the material will gradually resoften. Dislocations 
can be formed as single dislocations at the 
surface of a crystal or they can be formed as a 
dislocation pair inside the crystal. According to 
the calculations given in the paper about half as 
much energy is required to produce a certain 
amount of slip with surface dislocations as is 
required with interior dislocation pairs. Koehler 
has also calculated the strain energy stored in a 
highly work-hardened material. He assumed that 
the material then contains many dislocations 
arranged in a periodic two-dimensional array. 
The results for aluminum, copper, nickel, and 
iron are in good agreement with the experimental 
values. 

1 E. Orowan, Zeits. f. Physik 89, 634 (1934); M. Polanyi, 
Zeits. f. Physik 89, 660 (1934); G. I. Taylor, Proc. Roy. 
Soc. 145, 362 (1934). 

2F. Seitz and T. A. Read, J. App. Phys. 12, 100, 170, 
470, 538 (1941); W. F. Brown, Jr., Phys. Rev. 60, 139 
(1941); N. F. Mott and F. Nabarro, Proc. Phys. Soc. 52, 


86 (1940); F. Nabarro, Proc. Roy. Soc. 175, 519 (1940). 
3]. S. Koehler, Phys. Rev. 60, 397 (1941). 


New Books 


Research—A National Resource. II.—Industrial 
Research 











Report of the National Research Council to the 
National Resources Planning Board, December, 1940. 

Pp. 370+xi, Figs. 104, 2329 cm. U. S. Government 
Printing Office, Washington, D. C., 1941. Price $1.00. 
This ponderous report on “Industrial Research,” com- 
prising upwards of a quarter-million words and a hundred 
illustrations, printed on 370 large pages of a paper-bound 
book, and divided into seven major sections and numerous 
subsections, was submitted in compliance with a request 
of the National Resources Planning Board by the National 
Research Council. It was prepared under the supervision 
of a special survey committee of 26 members, of which 
F. W. Willard, President of Nassau Smelting and Refining 
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Committee, was chairman. The survey of industry was 
made by a specially appointed staff under the direction of 
Raymond Stevens, Vice President, Arthur D. Little, Inc. 
Each subsection was written by a different author. 

The authors define “‘Industrial Research” as ‘‘an organ- 
ized and systematic search for new scientific facts and 
principles which may be applicable to the creation of new 
wealth.” They estimate that in the United States 2200 
industrial research laboratories now employ over 70,000 
people, half of whom are professionally trained and half 
of the remainder technically trained. This is 41 percent 
more than were employed 2 years ago. Using the revised 
cost figure of $4000 per man-year, it is estimated that 
American industry is spending over $300,000,000 annually 
on research. 

The report discusses the purpose, scope and nature of 
‘industrial research,”’ its place in the industrial organiza- 
tion and in the national economy, cites successful examples 
of research in industry, presents charts and diagrams to 
show the location and extent of industrial research activity 
in the United States, compares American effort with that 
in the major foreign countries, and elaborates on the 
personal traits and disciplinary trainings which qualify 
men for careers in this field. 

The summary of findings in the report is followed by a 
set of recommendations: first, to industry, that it initiate, 
develop or expand its research program making use of the 
counsel and cooperation outlined in the report, and 
generously publish research results; second, to labor and 
industry, that they cooperate in studying fatigue and 
related matters of mutual concern; and finally, to govern- 
ment, that it foster industrial research as a national re- 
source contributing to the public welfare, support the 
publication of adequate abstracts of scientific and tech- 
nical literature, provide for revision and extension of 
critical tables of numerical data, and extend its research on 
standards of measurement. 

Noticeably missing from the report are: (1) a discussion 
of patents, which have been a powerful motivating in- 
fluence in industrial research; and (2) an explanation of 
the manner in which universities can and are engaging in 
practical and profitable research in cooperation with 
private research organizations. Although the authors 
acknowledge these omissions, they hardly justify them. 
Certainly, the important role being played by universities 
in defense research should warrant their more important 
mention as a national resource. Failure to give more 
‘consideration in the report to the possibilities of coopera- 
tive research in universities is hard to understand, since 
the authors attribute “the remarkable industrial growth 
which Germany experienced up to 3 or 4 years ago” to 
“the system of research in the universities and its coordina- 

tion with industry.” 

There can be little doubt as to the value of this govern- 
ment publication on “Industrial Research.” It is an 
encyclopedia of information, reviewing achievements, 
cataloging available facilities, and forecasting trends in 
this important field. The report should prove of value, not 
alone to the National Resources Planning Board and other 


770 


government agencies, but to any organization seekin- 
guidance on research matters or to individuals considerin, 
industrial research as a career. 

A. R. OLPIN 

Ohio State University Research Foundation 


Musical Acoustics 


By Charles A. Culver. Pp. 194+-xiii, Figs. 128, 1522 
cm, The Blakiston Company, Philadelphia, 1941. 
Price $2.50. 

This well-written and abundantly illustrated little work 
is the only American book in its field published within the 
last twenty-five years. Its style and level of presentation 
are carefully adapted to the needs of music students. 
All the essential topics are covered. The preface contains . 
particularly able statement of the reasons why musicians 
should study acoustics. A useful tabulation of symbols and 
terms follows. This reviewer wishes that the following 
definition had been included: 

Note: a printed or written character designating a musical sound of 
definite relative pitch; often misused in other senses. 

Most of the important modern findings are adequately 
presented, though the authoritative work of Stevens and 
Davis (Hearing, its Physiology and Psychology), Wiley, 
1938) seems not to have been consulted. Many original 
photographs of wave forms are shown. Some of these 
would be greatly improved if taken with a faster sweep so 
that no more than six or eight waves appeared on the 
oscilloscope screen. The system of relative pitch notation 
used by the author designates the lowest C on the piano 
as Co. A better choice would be C,, for the reasons set 
forth by Robert W. Young (J. Acous. Soc. Am. 11, 134 
(July, 1939).) 

It is unfortunate that the author’s clear statement of 
the nature of harmonic and inharmonic partials (p. 58) is 
not accompanied by a figure showing the harmonic series 
in musical notation up to perhaps the twentieth term. No 
mention is made of Doppler’s principle, nor is there any 
discussion of the nature and importance of transients in 
percussive and consonantal sounds. An appraisal of the 
requirements that must be satisfied by a really high fidelity 
loudspeaking system might well be given. 

The book is somewhat marred by the presence of minor 
errors, some twenty-three being noted by this reviewer. 
These include misspellings (mostly in words of foreign 
origin) and misstatements of fact. 

Despite its minor blemishes, the book so closely meets 
the needs of those who are teaching the subject that it 
merits general adoption as a classroom text. Any physicist 
who likes music but has not devoted much attention to 
acoustics will find fascinating reading here. One is tempted 
to append a similar promise for musicians who like physics; 
perhaps this is Utopian! Certainly any musician who wants 
to know more about the scientific basis of his art may turn 
to this book, assured that if the author cannot make 
acoustics easy, he can and will make acoustics plain. 

Cuas. WILLIAMSON 
Carnegie Institute of Technology 
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' Here and There 


Symposium on Non-Linear Circuit Theory 








Power and communication engineers are often con- 
fronted with problems involving non-linear circuit theory. 
The Basic Science Group of the American Institute of 
Electrical Engineers, New York Section, offers a Sym- 
posium of six lectures to acquaint electrical engineers with 
the material in this field; to present the more important 

, analytical and graphical methods for solving non-linear 
‘equations; and to stimulate further interest and research 


in this work. 


¥,< The lectures will be simplified and coordinated on the 


i 


assumption that the listeners have had no advanced 
mathematical training. The subjects, lecturers and dates 
are: 


; 1. General Nature of the Problem of Non-Linear Circuits and Out- 


: “line of the Various Methods of Attack. PRorrssor E. WEBER, Poly- 


= technic Institute of Brooklyn. Wednesday, November 5, 7 P.M. 
2. Power Series Solutions. Dr. L. A. MacCoit, Bell Telephone 
Laboratories. Wednesday, December 10, 7 P.M. 
3. Fourier Series Solutions. W. R. BENNETT, Bell Telephone Labora- 


- tories. Wednesday, January 7, 7 P.M. 


4. Graphical Methods. A. PrREISMAN, RCA Institute. Wednesday, 
. February 4, 7 P.M. 
5. Numerical Step-by-Step Method. ProressoR ALAN HAZELTINE, 
Stevens Institute of Technology. Wednesday, March 4, 7 P.M. 
6. Analytical Methods. Dr. P. Le CorBei__er, Lecturer, Graduate 
School of Engineering, Harvard University, and Lecturer, Electrical 
Engineering, New York University. 


Each lecture will be approximately one hour in length 
followed by a period of discussion. The fee for the complete 


* Symposium will be $1.50 for members of the American 


Institute of Electrical Engineers, American Society of 


'* Mechanical Engineers, American Society of Civil Engi- 


neers, American Institute of Mechanical Engineers, Insti- 
tute of Radio Engineers, American Mathematical Society, 
and American Physical Society, and $2.50 for non-mem- 
bers. Further information may be obtained from H. E. 
Farrer, American Institute of Electrical Engineers, 33 
West 39 Street, New York City. 


* 


Physicists Needed 


The number of junior physicists on the employment 
registers of the Civil Service Commission is not sufficient to 
meet the needs. A new examination announcement has 
therefore been issued under which applications will be 
accepted until further notice for junior physicist positions 
paying $2000 a year. The new examination is of the 
“unassembled” type, that is, applicants will not have to 
take a written test but will be rated on their education and 
experience. Persons who have received eligible ratings on 
the physics option of the 1941 Junior Professional Assistant 
examination need not file a new application as their 
eligibility will be continued on the new register. 

Junior physicists assist in or do research or investigative 
work in some field of physics. They are particularly needed 
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to work in the fields of radio and sound. Applicants must 
have completed a 4-year college course including, or 
supplemented by, at least 24 semester hours’ study in 
physics. In the rating of applications, extra credit will be 
given for additional undergraduate or graduate study in 
physics and for appropriate experience. Applicants must 
not have passed their thirty-fifth birthday. 

The Civil Service Commission also announces open com- 
petitive examinations for the positions of Senior Technical 
and Scientific Aid, $2000 a year, optional branches in 
radio and explosives; Technical and Scientific Aid, $1800 
a year, Assistant Technical and Scientific Aid, $1620 a 
year, Junior Technical and Scientific Aid, $1440 a year, 
optional branches in radio, explosives, chemistry, physics, 
metallurgy, and fuels. Duties will be to perform sub- 
professional technical and scientific work in the optional 
branches indicated; to conduct elementary chemical or 
physical testing or analysis; to make appropriate calcula- 
tions; and to perform related work as assigned. 

For some time the Commission has been accepting 
applications for physicist positions paying from $2600 to 
$5600 a year. Physicists interested in securing Government 
employment are urged to seek further information about 
these examinations from the Commission’s representative 
at any first- or second-class post office or from the central 


office in Washington, D. C. 
* 


Necrology 


Dr. Harry S. Hower, Professor and Head of the Depart- 
ment of Physics at Carnegie Institute of Technology, 
Pittsburgh, Pennsylvania, fellow of the American Physical 
Society and of the American Association for the Advance- 
ment of Science, died on October 10 at the age of 64. 

Dr. Samuel Imes, Professor and Head of the Department 
of Physics at Fisk University, Nashville, Tennessee, died 
on September 12 at the age of 58 years. 

Mr. C. C. Potwin, acoustic consultant for Electrical 
Research Products, Inc., New York City, Fellow and 
Treasurer of the Acoustical Society of America, died on 
September 25 at the age of 33. 

Dr. Carroll Mason Sparrow, Professor of Physics at the 
University of Virginia, died on August 30 at the age of 61. 


* 


Buhl Foundation Grant 


A fund of $40,000 per year over a 5-year period has been 
made available to the Departments of Biology, Chemistry 
and Physics at the University of Pittsburgh, primarily for 
research fellowships and supplies, according to a recent 
issue of Science. Three-fourths of the fund is contributed by 
the Buhl Foundation and one-fourth by the University. In 
the division of physical chemistry, specific heats, heats of 
combustion and heats of dilution will be studied in relation 
to the molecular structure of sugars, fats and amino acids. 
Research involving the use of radioactive tracer elements 
and spectroscopy will be carried out in cooperation with 
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the Department of Physics, where additional work will 
be supported in the field of atomic physics (a cyclotron is 
now under construction) and in the use of the electron 
microscope. 


* 


Recent Appointment 


Edwin H. Colpitts has been appointed Director of the 
Engineering Foundation effective October 1, 1941, succeed- 
ing the late Dr. Otis E. Hovey, according to an announce- 
ment of the Foundation. Established in 1914 by a gift of 
the late Ambrose Swasey, the Foundation is a joint cor- 
porate agency of the four major engineering societies, 
namely, the American Society of Civil Engineers, the 
American Institute of Mining and Metallurgy, the Ameri- 
can Society of Mechanical Engineers and the American 
Institute of Electrical Engineers. The Foundation has no 
laboratories or technical staff, but sponsors and supports 
investigations in university, government and _ private 
laboratories. Dr. Colpitts, who is an LL.D. of Mount 
Allison University, is a Fellow in the American Physical 
Society, in the Acoustical Society of America, in the 
American Institute of Electrical Engineers, and in the 
Institute of Radio Engineers. From 1933 to 1937 he was 
Vice President of Bell Telephone Laboratories. 


* 
Physicists in Defense 


Professors W. Warren Stifler and Theodore Soller have 
been given leave of absence from the Department of 
Physics at Amherst College to enable them to engage in 
national work. Professor Soller will undertake 
research at Massachusetts Institute of Technology while 
Professor Stifler will conduct special courses in optics at 
Columbia University. 


defense 


* 
Program in Applied Mechanics 


The nation’s first full-time program of instruction and 
research in applied mechanics began its Fall session at 
Brown University on October 2. The courses offered deal 
with specialized problems in mathematics which bear upon 
aviation and airplane construction, with shipbuilding, gun 
manufacturing, armor-plating on tanks and armored cars, 
optical instruments, and other weapons and instruments of 
war requiring a high degree of technical knowledge to 
plan, design and make. 

Dean Richardson announced the following have been 
appointed to the teaching staff this year: John L. Synge, 
Professor of Applied Mathematics at the University of 
Toronto, who is engaged in research for the National 
Research Council of Canada; Richard von Mises, Head of 
the Department of Aeronautics in the Graduate School of 
Engineering at Harvard, formerly Director of the Institute 
for Applied Mathematics at the University of Berlin; 
Jacob D. Tamarkin, Professor of Mathematics at Brown 
and Editor of the international journal, Mathematical 
Reviews; Willy Feller, Assistant Professor of Mathematics 
at Brown and formerly of the Universities of Kiel, Copen- 
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hagen and Stockholm; Stefan Bergman, expert in applied 
mechanics, who was formerly an Instructor at the Institute 
for Applied Mathematics at the University of Berlin and 
at the Technological Institute of Tomsk; Willi Prager, 
Professor of Engineering at the University of Istanbul and 
formerly associated with the Aerodynamical Institute of 
Gottingen. 


* 
Bulletin of Mathematical Biophysics 


The December, 1941, issue of the Bulletin of Mathe 
matical Biophysics, edited by N. Rashevsky and published 
by the University of Chicago Press, will contain the 
following articles: 

A Neural Mechanism for Discrimination: IV. Monocular 
Depth Perception, HENRY STANTON 

Some Preliminary Considerations Concerning Concen- 
tration of Oxygen in Tissue, INGRAM BLOCH 

Note on a Possible Application of Some Concepts of 
Topology to Asymmetric Organization of Protoplasm, G. F. 
GAUSE 

Statistical Distribution of Impedance Elements in Bio- 
logical Systems, A. M. WEINBERG AND A. S. HOUSEHOLDER 

A Theory of Steady-State Activity in Nerve-Fiber Net- 
works: III. The Simple Circuit in Complete Activity, 
ALSTON S. HOUSEHOLDER 

Distribution of Response Times, H. D. LANDAHL 
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* 
Honors Conferred 


Honorary degrees have been conferred by the University 
of San Marcos at Lima and by the University of Arequipa 
upon Dr. Arthur H. Compton, Chairman of the Depart- 
ment of Physics at the University of Chicago, it is noted in 
Science. Dr. Compton was also made an honorary member 
of the Peruvian Chemical Society, and a corresponding 
member of the Brazilian Academy of Sciences and of the 
Academy of Exact Sciences at Lima, during his recent 
visit to South America as director of the Andean Cosmic- 
Ray Expedition. 

It is also noted that, at the convocation on September 
29, following the symposia on ‘‘New Frontiers in Education 
and Research” held at the University of Chicago on the 
celebration of its fiftieth anniversary, in cooperation with 
the American Association for the Advancement of Science, 
the following were among those receiving degrees: E. O. 
Lawrence, Professor of Physics at the University of Cali- 
fornia; Robert A. Millikan, Chairman of the Executive 
Council of the California Institute of Technology; Linus 
C. Pauling, Professor and Chairman of Department of 
Chemistry, California Institute of Technology; and Henry 
N. Russell, Director of the Princeton Astronomical Ob- 
servatory. 
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Matrices, Tensors, or Dyadics for Studying Electrical Networks? 


MyriL B. REED 
Illinois Institute of Technology, Chicago, Illinois 


(Received April 18, 1941) 


Dyadic algebra, because it has been designed for. three-dimensional systems, must undergo 
further development to be effective for more than three-mesh networks. There is little point in 
carrying out such a development, since the much used cross product and the symbolism of 
dyadics both lose their simplicity for more than three dimensions. Of the outstanding properties 
of tensor analysis, only the special technique of linear transformations has been used in electrical 
network studies. But linear transformations are basic to the theory and manipulations of matrix 
algebra. Therefore, the simplest and most effective approach to network transformations is 
through matrix algebra. The primitive-network impedance matrix of tensor treatments is in 
reality nothing more than the matrix of the branch impedances of a network, and the linear 
transformation employed to formulate the mesh-current impedance matrix is simply the trans- 
formation which relates mesh and branch currents. Furthermore, the “interconnection of net- 
works”’ transformation is merely a transformation from more than the minimum number of 
mesh currents to the minimum number. Matrix transformations may be used to establish 





networks which are equivalent. 


INTRODUCTION 


YADICS,! tensors,2* and matrices**® have 

all been expounded at considerable length, 
with the claimed virtues and powerfulness of 
these methods forming a prominent part of the 
discussions. The claims, of almost all-embracing 
power of analysis, particularly of tensors, have 
not gone unchallenged.*’ The present paper has 
been written in an. endeavor to show how far 
these claims for unusual merit are justified, and 
to show that the matrix method is the most 
useful. 

DYADICS 


The method of dyadic algebra was developed 
by mathematicians to facilitate the study of 





1A. P. T. Sah, Dyadic Circuit Analysis (International 
Textbook Company, 1939). 

2G. Kron, Tensor Analysis of Electrical Networks (John 
Wiley and Sons, 1939). 

3L. V. Bewley, “Tensor algebra in transformer circuits,” 
Elect. Eng. 55, 1214 (1936). 

4*L. A. Pipes, ‘‘Matrices in engineering,’’ Elect. Eng. 56, 
1177 (1937). 

5 M. B. Reed, ‘‘Properties of three-phase circuits deduced 
with the aid of matrices,’’ Elec. Eng. 57, 74 (1938). 

6 J. Slepian, ‘“‘A parable on tensor analysis,” Electric J. 
33, 541 (1936). 

7W. H. Ingram, “‘Tensors and electrical networks,”’ 
Elect. Eng. 58, 274 (1939). 
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physical systems which are most’ effectively 
represented in three-dimensional space. This fact 
suggests that dyadic algebra will, therefore, be 
restricted in its applicability to three-dimensional 
forms or, in electric circuit parlance, to three- 
mesh networks. Three-phase systems are assumed 
to be particularly receptive to dyadic algebra 
treatment, but, in actual fact, not all three-phase 
systems can be handled as three-mesh networks, 
as for example, a delta-delta system—a four- 
mesh network. 

The extension of dyadic algebra to more than 
three dimensions is hardly justifiable. The 
symbolism is more complex than matrix sym- 
bolism and the cross product, which holds a very 
prominent position in dyadic algebra, loses its 
simple determinantal form and becomes very 
complex beyond the third order. As a result, the 
dyadic method has been generally abandoned by 
mathematicians in favor of matrix algebra—a 
simpler symbolism with no restrictions as to 
dimensionality. 


TENSORS 


The method of tensor analysis was developed 
by mathematicians in an effort to derive mathe- 
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matical forms which are the same for all coor- 
dinate systems. The outstanding value of tensor 
analysis, as conceived by these mathematicians, 
is the use which can be made of the fact that if 
the components of a tensor are zero in any coor- 
dinate system, they are zero in all coordinate 
systems. This outstanding feature of tensor 
analysis has not been used in electrical network 
studies. Of the several different features of tensor 
analysis, the most widely known are the index 
notation and the technique of transforming 
relations from one coordinate system to another. 
It is these two features for which so much is 
claimed in the solution of problems in electrical 
network analysis. 

The transformation of liaear systems from one 
coordinate system to another, while an integral 
part of tensor analysis, is also fundamental to 
the theory of linear equations and matrix algebra. 
This fact will be discussed later. The index nota- 
tion has not been of particular use in electrical 
network studies because the formulations of 
tensor analysis in which the index notation is 
particularly helpful do not appear in network 
studies. For example, the formulation of the 
Christoffel symbols is greatly facilitated by the 
index notation but these symbols have not 
appeared useful in network analysis or synthesis. 
The distinct advantage of index notation over 
other forms in the formulation of m-way tensors 
(n>2) disappears in the manipulation of the 1- 
and 2-wavy forms of electrical networks. In fact, 
the advantage is now in favor of matrix sym- 
bolism. 

One other point in connection with tensor 
analysis should be mentioned. It is claimed that, 
from a so-called primitive network of n branches, 
any n-branch network can be established by a 
linear transformation. It will be shown in what 
follows that this statement of the situation is 
incorrect and an interpretation, expressed in 
matrix form, will be given to the assumed tensor 
transformation from a primitive network. 


MATRICES 


Fundamentally, all electric network problems 
must be formulated from five basic laws—the 
laws relating current and voltage in and across 
a resistor (Ohm’s law), across an inductor (Lenz’s 
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or Faraday’s law), and across a capacitor, and 
Kirchhoff’s voltage and current laws. The use of 
these laws leads to the system of simultaneous 
equations, which are ordinarily idealized into the 
linear equations used for circuit analysis. Basi- 
cally, therefore, the properties of electrical 
networks must all be deduced from a system of 
equations which are most simply and easily 
symbolized by matrices.**!° The completely 
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extraneous concepts of phases as space dimen- 
sions, primitive networks, covariant and contra- 
variant variables, index notation, and n-way 
tensors never arise. Instead, the properties of an 
electric circuit follow naturally and simply from 
the properties of matrices which have been 
developed by mathematicians for the express 
purpose of establishing the properties of systems 
of linear equations. One curious fact is evident 
in the writings of tensors and dyadics; that is 
that matrices and the properties of matrices are 
widely used to develop the results which are 


.supposed to arise from tensors and dyadics."? 


PRIMITIVE NETWORKS 


In line with the concepts of tensor analysis, 
the tensor treatment of the electrical network 
has been formulated on the basis of transforming 
from the coordinate system (primitive network) 
which permits the simplest expression for the 


S$L. A. Pipes, “Transformation theory of networks,” 
Elect. Eng. 59, 123 (1940). 

*L. A. Pipes, “The operational calculus. I1,”’ J. App. 
Phys. 10, 258 (1939). 

10 L. A. Pipes, “Matrix theory of oscillatory networks,” 
J. App. Phys. 10, 849 (1939). 
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current and voltage relations to another coor- 
dinate system in which the network under con- 
sideration may be expressed. However, there are 
some outstanding difficulties in connection with 
such an interpretation. 

The network of Fig. 1(b) will be considered. 
The so-called primitive network of this circuit is 
given in Fig. 1(a). This primitive network is 
assumed to be the basic network from which all 
six-element networks can be deduced by linear 
transformations of a general form, of which the 
following is a particular example for the networks 
of Fig. 1 


lz] 1 0 of 
I3/\| |\1 0 -1 


\Zs\) 
arn Ze _O 8 ANN ane 
\Ip|| = id _ 0 me 0 Is = iC! lI ml. (1) 
I; 0 —-1 1\| '"#!! 


Ie] 0 O —1 


Similar transformations can be established for 
other six-element networks. Such transforma- 
tions, however, must have the restriction, ignored 
by tensor analysis proponents, that all networks 
deduced from any particular primitive network 
must have the same mutual induction between 
the same impedances as the primitive network. 
Actually, therefore, the so-called primitive net- 
work could be primitive only in a restricted sense. 

Furthermore, even a superficial examination 
of the so-called primitive network and the 
network derived therefrom shows that power 
absorption in the two networks is not the same, 
so that the power invariance under the trans- 
formation of Eq. (1), assumed in the tensor 
treatment, loses its meaning. Thus, the trans- 
formations deduced by means of Eq. (1) under 
the assumption of power invariance cannot be 
justified. Actually, the transformation of Eq. (1) 
is valid if based on the proper premises and does 
lead to valid conclusions, as will now be shown. 

If a system of equations is written which gives 
the branch currents of a network in terms of the 
mesh currents of that same network, the resultant 
matrix equation will have the form given in Eq. 
(1). In fact, Eq. (1) is the system of equations 
which expresses the branch currents of Fig. 1(b) 
in terms of the mesh currents of that same figure. 
Assuming a different set of mesh currents will, 
of course, require a different c matrix. Using 
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instantaneous forms for currents and real number 
c matrices, the branch currents and mesh currents 
for any network will be related by 


| Zo(¢) || = || Cl] || Zm(2)||. (2) 


The voltage and impedance transformations 
may be established by making use of the obvious 
fact that the power delivered by the generators 
is equal to that absorbed by the network. Thus, 
the instantaneous power delivered to a network 
may be expressed as 


p(t) = || Em(t)|| || Zm(t) ||, (3) 
where, for the network of Fig. 1(b), 


|ex(t) —es(¢)|| ||e1(¢)| 
p(t)=|| —e2(t) I \ia(t)| 
e3(t) I e\lis(t) 
i1(t) 
= |/e1(t) —e3(t), —e2(t), es(t)|||\72(2)||. (4) 
73(t)|| 


The power delivered by the generators may also 
be expressed as the sum of the products of the 
generator voltages and the corresponding branch 
currents; i.e., 


p(t) = ||Eo(t)||!|Lo()|!, (5) 
which, for Fig. 1(b), is 





| ex(t)||| | int), 

| — €3( (t)||, iso) 
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Finally, the power absorbed by the network may 
be expressed, in terms of the branch currents and 
voltages, as simply the sum of all the corre- 
sponding branch current and voltage products; 
thus 


p(t) = || Vo(t)|| | Zo(2) |], (7) 
which, for the network of Fig. 1(b), is 








Ilv (| |é1' (t) 
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The instantaneous power expressions of Eqs. 
(3), (5), and (7) may be equated, even though 
the current and voltage matrices are themselves 
not equal, because they are simply different 
expressions for the power absorbed by the 
network. Hence 


|| Em (t)!! || Ln (t) |! = |! Eo (t) |! «| Do (t) 
=|) Vo(t)! | D,(2) ||. (9) 
The voltages may be related from these equalities 
and Eq. (2). 
If the two expressions on the left of Eq. (9) 
are considered and if Eq. (2) is used, 
|| Eon (£) | |] Lm (¢)!| = || Eo (¢) |{ ll Cll | Zm (2) |, 
so that 
[|| Bm (t) || e— || Be(t) |! el] Cl] JI Zm(¢) || =0. (10) 


Then, since the matrices in the brackets are 
independent of |\I,,(t)||, the bracketed term must 
be zero since this equation is valid for all values 
of time. Thus 


|| Em (t) || e= || Eo (t)!| | Cl (11) 

and since, in general, 
(!|A || \|Bl| Je= |] Ble A ll, (12) 
|| En (t)!! = || Cl] «|| B,(6)!!. (13) 


This equation is also valid for the complex 


Z1+Z3 0 
LSet Ea =e 


This matrix is, of course, the usual mesh current 
matrix which arises by writing the mesh voltage 
equations directly from the circuit. 

The transformations just worked out show 
that, by using only the simplest concepts of 
matrix algebra, conclusions can be reached by 
simple and logical applications of the basic laws 
of electric circuits; whereas, even with all the 
“excess baggage’’ of tensor analysis, these same 
conclusions have not been reached logically. 

The formulas developed in the foregoing may 
have value for computational purposes.’ The 
impedance formula of Eq. (19) is probably the 
most useful, since writing ||Z,|| from the circuit, 
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Z2+Z5 


so that 
|Em|| = || Cl] || Eel]. (14) 


In an identical manner, from the first and last 
terms of Eq. (9), it follows that 


E,,|| = Cll «| Ve . (15) 


Turning now to the relation of the currents and 
voltages by means of the impedances, the branch 
voltage drops are 


| Voli =|! Zell || Dell, (16) 


which, for Fig. 1(b), are (complex number form, 
with mutual induction between Z; and Z,) 


Vy'|| |Z, 0 O O O OZ,’ 
| 10 Zs Ze O O ON IITs’ 
Vi)_|0 Zy Z, 0 O O I l,'1| (17) 
Ve) 10 O O Ze O OfliiZe'\\ 
V1 10 0 0 O Ze ONT,’ 
Vell 10 0 0 OO O Zell |iIe 


If Eq. (16) is multiplied on the left by ||C|!, and 
if ||J,|| is replaced by its equal from Eq. (2), 
|| Em || = || C'l «|| Vol] =] Cll el Zo) ||C | ||Zm|], (18) 
from which, as a matter of convenience, we may 
write 
Zm|l = || Cl el]Ze|| || C (19) 


For the network of Fig. 1(b), this impedance 
matrix is 


—Z3t+Za3 | 


| 
—Zs . (20) 


Z3t+ZitZst+Ze—2Z,3) 


and at the same time keeping the mutual- 
induction factors in their proper positions, is 
simple; whereas, in the final form, as ||Z,,||, the 
mutual-induction coefficients are very easily 
omitted or incorrectly placed. The branch 
currents and voltages are easily obtained from 
Eqs. (2) and (16), or merely by inspection of the 
network, as ||C!! was determined in the first place. 





INTERCONNECTION OF NETWORKS 


The concept of tensor analysis has been used 
to form a transformation, under power invariance, 
which is said to express the behavior of isolated 
networks when they are connected together. For 
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example, the behavior of the network of Fig. 1(b) 
































or Fig. 2(c) is expressed in terms of the two 
isolated networks of Fig. 2(a) and Fig. 2(b). 
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But the assumption that the power absorbed 
by the networks of Figs. 2(a) and (b) is the same 
as that absorbed by the networks of Figs. 1(b) 
or 2(c) is obviously not generally correct. The 
formulas must, therefore, be derived from a dif- 
ferent assumption to be valid. As will now be 
shown, the transformation for “network inter- 
connection”’ is really a transformation from the 
branch aspect of a network to the mesh aspect, 
where more than the minimum number of mesh 
currents are employed, followed by a trans- 
formation which reduces the mesh equations to 
the minimum numer. The network of Fig. 2(c) 
will be used to illustrate these transformations. 

First, the transformation from the branch 
aspect to extra-mesh- current aspect is—Eq. (2) 





| 
" i" | non | | ||| Cir’ t 
Zn || = ty t Z\\\|C’|| = ; 


0 1Coo' I ,i{|| 0 


—using compound matrices to make the sig- 
nificance of the transformation more evident, 
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“| 0 ll Coo’! I! Tao | ax | | C"| 1? (21) 


yy 


Hil 


Second, the transformation from the extra- 
mesh-current aspect to the minimum-mesh- 
current aspect is 











11 O Oll,_, 
Z| 
g , 10 O 1) {07I 
| Zn’ || a ! ha \| Zon! =] =|} || Zell 
JO 1 lll | 
|0 0 i ™ 

“li “ie |Znll. (22) 

Finally, the branch currents, in terms of the 


minimum number of mesh currents, is, from Eqs. 
(21) and (22), 


| Lol] =|] ||| Lm’|| = 1] C7] || C77} || Z| 


This combined transformation 
transformation given by Eq. (1) 

In order to show how the foregoing trans- 
formations might possibly facilitate computa- 
tions, Eq. (19) will be expressed in terms of the 
transformation of Eq. (21). Thus, if the branch 
aspect impedance matrix is expressed in com- 
pound form 


(23) 


is iaiatile the 


HII Z va| 0 


Zeal\|ll| 0 (24) 


It may be possible to further subdivide the matrices in any particular case. In fact, the subdivision 
may be carried to any extent desired as long as the diagonal form given in Eq. (24) results. It should 
be noted that there can be no mutual induction between any of the network elements contained in 
||Z11|| and those contained in ||Z22!| of Eq. (24). 

Performing the multiplications indicated in Eq. (24) gives 


|Cir'|| 0 


| |Cus"[lellZaa (25) 
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The meaning of these resultant matrices can best be shown by a particular example. For the network 
of Fig. 2(b), Eq. (24) is 


'Z, 0 0O]0 O Oo] 
) 0 Zs; Z30 0 0 








az.’ = | Curl 1] / | : ott. nd : : Hea | . po (26) 

| 0 Coo" |\¢ 0 O 0 » 0 O 0 | Coe 1} 

, 0 0 0 0 Zz o|l | 

10 0 0 |0 O Ze 
Substituting the c matrices from Eq. (21) into Eq. (26) and multiplying gives 

| ZitZs, —ZstZas | P | 
ZH al —Z3+Zu3 Z3+Z4—2Z,3| - — |Z 11" || 0 | (27) 

aris | 0 \Zot+Zs —Zs 0 |Z 20" Ae 

| —Zs Z5+Ze 

The submatrices of this expression are the usual The foregoing discussion has shown the basis 


mesh current matrices of the circuits of Figs. on which a transformation of currents may be 
2(a) and 2(b). But the combined behavior of performed and has also shown that intercon- 
these two networks in no way corresponds to nection of networks cannot be performed as the 
the behavior of the network of Figs. 1(b) or tensor treatment has indicated. 

2(c). Certainly the power absorbed is not the 

same. Therefore, the transformation is not one TRANSFORMATION OF ONE NETWORK TO 

of power invariance. The question thus arises as ANOTHER EQUIVALENT ONE 

. to just what the transformation of Eq. (21) does 
signify. The answer is simple. First, there is no 
doubt that Eq. (17) represents the branch aspect 
of any network having six impedances, however 
arranged, with the ones designated as Z; and Z, 
mutually coupled. The transformation of Eq. 
(21), therefore, is the transformation from the 
branch aspect of Figs. 2(a) and 2(b) to the mesh 
current aspect of these same two networks, com- 
bined into a single compound matrix consisting 
of independent submatrices. On the other hand, 
the transformation of Eq. (1) is a transformation 
from the branch currents of Fig. 1(b) to the ; Sraee 
ine emanate ue wet endl tx tee formation and assuming it does represent the 
snene tibial of temdionmetion os cither ef the SON a te two sets of mesh currents, because 
submatrix transformations of Eq. (21). Trans- /« and J, of Figs. 2(a) and 2(b) are not neces- 
formations of the kind used here must, therefore, sarily equal to I; of Fig. 2(c). With this error at 
be interpreted as a transformation from branch the beginning, of course, the remainder of the 


As has already been shown, the transformation 
given by Eq. (22) is actually a transformation 
under invariance of power only if the currents 
all belong to the same network, as in Fig. 2(c). 
In this circuit the currents J;, J,, and J, are 
obviously all equal and the transformation, even 
though singular, keeps power invariant. How- 
ever, if this transformation of Eq. (22) is viewed 
as the transformation from the mesh currents 
of Figs. 2(a) and 2(b) combined, to the mesh 
currents of Fig. 2(c), it is mot correct; i.e., an 
error is made initially in setting up the trans- 


to mesh aspect of networks. transformations and any resultant conclusions 
The transformation of Eq. (1), when expressed re incorrect. 
as the two transformations of Eq. (23), mav It 1S possible to establish networks which are 


possibly be of use in reducing the size of the equivalent to a given network by transformations 
matrices to be handled. For example, Eq. (27), similar to the foregoing. The matrices should, in 
or its equivalent, was formed by handling two general, be real and non-singular. Then, for any 
second-order matrices in place of one third order. given n-mesh network, ‘on assuming any such 
As the number of meshes increases, whatever non-singular real matrix and using Eqs. (1), (14), 
advantage this scheme may have will increase. and (19) (where b is changed to some other 
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appropriate symbol, say m,) another n-mesh 
impedance matrix may be established and the 
circuit, or circuits, which this derived impedance 
matrix represents will be equivalent to the 
original network insofar as the instantaneous 
power delivered by the generators is concerned. 
The system of generators in the two equivalent 
networks will, in general, be different. 


CONCLUSIONS 


The foregoing discussion has shown that the 
concepts of tensor analysis are not necessary for 


the complete understanding and development of 
the transformations which have been stressed as 
arising solely from the tensor concept. Further; 
it has been shown that the tensor concept with 
its primitive network does not correspond to 
physical reality and that the transformation 
used can only be justified from an entirely dif- 
ferent viewpoint. This viewpoint is supplied 
easily and naturally by matrix methods because 
such methods have been developed to facilitate 
the manipulation of linear equations, which are 
basic in electrical network studies. 





The Evaporation of Molten Metals from Hot Filaments 


WaLLace C. CALDWELL 
Physics Department, Iowa State College, Ames, Iowa 


(Received May 12, 1941) 


The selection of filament material from which another metal can be evaporated successfully 
depends largely upon the three properties of wetting (W), reaction or alloying (R), and evapora- 
tion (EZ). Nine filament materials and twenty-seven evaporated metals have been investigated 
and the results summarized in a table using the above symbols with subscripts 1, 2 and 3 to 
indicate the degree of favorable properties. A second table lists satisfactory filament materials 


for the evaporation of twenty-six metals. 


EVERAL investigators! have used the evapo- 
ration technique for making surface mirrors, 
thermocouples, high resistances, and other de- 
vices. The method usually used involves attach- 
ing the metal to be evaporated to a refractory 
filament that is heated electrically. Tungsten 
was nearly always used as the filament until 
Countryman? found that for some metals other 
filament materials were to be preferred. In the 
present investigation a determination of the best 
filament to be used with each metal was made. 
The apparatus consisted of a bell jar ground to 
fit a flat glass plate and the associated mercury 
diffusion pumping system. Necessary electrical 
leads and supports were sealed in the base. 
A mercury pool sealed the joint between bell jar 
and plate. 


1C. Hawley Cartwright, ‘‘Notes on the evaporation of 
Ag, Be, Cr, and Si,” Rev. Sci. Inst. 3, 298 (1932); C. 
Hawley Cartwright, and John Strong, “An apparatus for 
the evaporation of various materials in high vacua,” Rev. 
Sci. Inst. 2, 189 (1931); John Strong, Procedures in Experi- 
— Faas (Prentice-Hall, New York, 1938), pp. 
168-187. 

2M. Alden Countryman, ‘‘Notes on the wetting of 
filaments by molten metals,” J. App. Phys. 8, 432 (1937). 
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If the metal to be evaporated was available in 
the form of wire, it was wound tightly about the 
filament at the bottom of a couple of sinuous 
loops; if it was available only in lump form, the 
lumps were placed in a vertical conical helix. 

In order that a metal evaporate from a fila- 
ment satisfactorily the metal must adhere to the 
filament, must have sufficient vapor pressure to 
evaporate while still at a temperature below the 
melting point of the filament material, and must 
not form an alloy with the filament which has a 
melting point below the temperature at which 
evaporation will ensue. A knowledge of the 
degree to which the filament material meets 
these requirements makes possible the selection 
of the filaments most suitable for the evaporation 
of each of the metals investigated. 

In Table I are given the results of visual 
observations of wetting of the filament by the 
molten metal, of evaporation, and of any obvious 
alloying of the molten metal with the filament 
material. An explanation of the abbreviations is 
given in the legend accompanying the table. 
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Tungsten W, W; W, W, We W; M; M; Ws; We W, W; 
3370-—5900°C mn 1 E, E, EB, es E; E; t3 Ey 3 2 
Re R, R, R Ri Ri R; R; Ri Ri R2 l 
Tantalum W, W. W, W, W. W; Ws; M; W:2 W, W, M; W; 
2850 4100°C on 1 Oe on on Ey 1 E; E; E, Ey E, Ey; 1 
R» R, R, R, R, R, R; R, R, Ri R2 R; Ri 
Molybdenum W, Ws; W, W, W.2 Ws; W2 M; W2 W, W.2 Wi W; 
2620-—3700°C on es E, | Oe on E; E; E; Ey Ey a1 E; 1 
R» R, R, R, Ri R, R; R, R, R,; R, R; Ry 
Columbium W, W; W, M; W2 W, W.2 2 W,2 W, Wi W; 
1950-—2900°C or 3 oF E, om Ey E, EF, E2 Es; 1 
R: R, R, R, R, R, R; Rk R2 R2 R; 
Platinum W, Wi 
1773-4300°C E; v3 
R; R; 
Iron W, We W, M; Ws; W2 W, W, W, 
1535 3000°C Es; On om E. l Ey Ey 3 Ey 
R; R; R, Ri Ri R, R; R> R, 
Nickel Wi W, W, W, W, Ws; W, W, W;, 
1452-—2900°C Ey v3 FE, E; E3 E; E, 23 Ei 
Re R, R, R; R, R, R3 Re Ri 
Chromel W, W, W, M; W, W, W, W, Wi 
1350°C Es; Ey E; E; E, Ey Ey 3 E, 
R; R; R, R, R, R, R; R2 Ri 
Silver W, W2 W, We W; W, 
960-—1950°C E; E; E, E; E: E; 
R; Rs Ri R3 Ri R; 


* W:—The metal wets the filament well. 
W:—The forces of wetting are not large. Drops of the metal may occasionally fall from the filament. 
W:s—The metal will not adhere to the filament. 

M;s—The metal is not melted. 
E,— Evaporation of the metal takes place readily. 
E2—Evaporation of the metal proceeds slowly. 
Es—No evaporation of the metal takes place. 
R:—No reaction of the metal with the filament. 
R2—Reaction of the metal with the filament proceeds slowly usua!ly without much effect on evaporation. 
R:—Reaction of the metal with the filament results in burning out the filament. 
B—The filament burns out due to heating the metal to nearly the melting point of the filament. 


NOTES 


1. In some cases where the metal does not wet the filament, the surface tension is such that the globule of molten metal is supported mechanically 
by the spiral filament. Although evaporation occurs, the occasional dropping of molten metal from the filament is an undesirable characteristic. 

. In an attempt to evaporate columbium, wire with a diameter of 0.03 in. was used. Apparently the heat transfer from filament to the columbium 
was not efficient enough to melt the columbium, but fine columbium wire wrapped tightly about the filament should melt quite easily. 

3. Iron is difficult to evaporate because it forms low melting point alloys. However, if several turns of tungsten are wrapped tightly about a 
short length of iron wire, the iron can be evaporated successfully. An iron-tungsten alloy is formed along the whole length of the helix, but a core 
of tungsten is left unaffected so that the filament is maintained intact. Differential evaporation of iron from the molten tungsten-iron alloy probably 
takes place. Nickel can be evaporated by this same procedure. 


4. Due to its low vapor pressure, platinum evaporates very slowly. If very fine platinum wire is wrapped tightly around a tungsten filament, the 
efficiency of heating is adequate to melt and to evaporate the platinum. 


5. A shell of some refractory compound of strontium separated the molten strontium from the filament so that the wetting properties could not 
be examined. 


Occasionally, when additional information should 
be given for a certain metal-filament combina- 
tion, a number is inserted in the table and the 
information given in the notes. Under the name 
of each metal to be evaporated at the top of 
the columns, and under the name of each filament 
listed in a column at the left, are given the 


melting and boiling points at atmospheric pres- 
sure. The metals to be evaporated are listed 
alphabetically, the filament materials in the order 
of decreasing melting points. 

Consultation of Table I enables one to select 
the best filaments for the evaporation of each 
of the metals considered in this investigation, 


780 JOURNAL OF APPLIED PHYSICS 








Tu 


33) 





.a™Maw \ 


TABLE I.—Continued. 
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Tungsten M: W W, Wi W3 W2 W, W, W2 W, Ws; W, W, W, 
3370-5900°C om Ey 3 4 1 Ey a1 Ey on E2 EB; Ey Ey EE 
Ri Ri R; R2 Ri Ri Ri Ri Ri Ri Ri R; R2 Ri 
Tantalum M; W, W, W, W3 W, 5 W, W2 M; W. WwW; Wi Wi 
2850-—4100°C om a1 3 E3; 1 Ey Ey on on E; a1 Ey D3 om 
R Ri R; R; Ri Ry R, R, R, R R, R, R; R; 
Molybdenum M; Wi W, W, VW 1 W, if! 1 W, W3 M; U W, W, U 1 
2620-3700°C Ey EF, o3 E; E, Ey E Ey 1 Ey Ey E, EB Ey 
R, R, R; R; Ry Ry Ry R, R, Ry R, R; R, R, 
Columbium M3; W, W, W, W, W, W, W, W, Ws; W, W, Wi 
1950—2900°C E, on 23 E3 Ey EB, Ey Ey EB, 1 E2 EE, E, 
R, Ri R; R; Ri Ry Ry R, R, R, R, R, R, 
Platinum W, 
1773-—4300°C E3 
R; 
Iron M; M; W, W, W, W, W, We W, 
1535-—3000°C Ey E, EF Ey Ey Ey EB Ey E; 
R, B Ri R, Ri R, R Ri R; 
Nickel M; W, W2 W, 5 W, W, Ws; WwW, 
1452-—2900°C Ey E; Ey Ey Ey Ey Ey E; D3 
Ry B R, Ri Ri R, Ri Ri R; 
Chromel M; Wi W, W, Wi W, W3 W, W, 
1350°C Ey E. Ey Ei Ey Ey 1 Ey E; 
Ri B R, R R, R R, Ri R; 
Silver W, M;3 W, W, Ws; W, 
960—1950°C Ey E3 E; E; 1 Es; 


TABLE II. Satisfactory filament materials for the evaporation of various metals (listed in order of merit). 


R: 








Ri Rs 








R3 Rs Re 





METAL FILAMENTS METAL FILAMENTS 
Aluminum Tungsten, tantalum, molybdenum, colum- | Manganese Tungsten, tantalum, molybdenum, colum- 
bium bium 
Antimony Chromel, tantalum, tungsten Nickel Tungsten 
Barium Tungsten, tantalum, molybdenum, colum- | Platinum Tungsten 
bium Selenium Chromel, iron, molybdenum 
Beryllium Tantalum, tungsten, molybdenum | Silver Tantalum, molybdenum, columbium, iron 
Bismuth Chromel, tantalum, tungsten | Strontium Tungsten, tantalum, molybdenum, colum- 
Cadmium Chromel, columbium, tantalum bium 
Cobalt Columbium | Tellurium Tungsten, tantalum, molybdenum, colum- 
Copper Columbium, molybdenum, tantalum bium 
Germanium Tantalum, molybdenum (other filaments | Thallium Nickel, iron, columbium, tantalum 
not tried) Thorium Molybdenum 
Gold Tungsten, molybdenum Tin Chromel 
Iron Tungsten Titanium Tungsten, tantalum 
Lead Iron, nickel, chromel Vanadium Tungsten, molybdenum 
Magnesium Tungsten, tantalum, molybdenum, colum- | Zinc Tungsten, tantalum, molybdenum, colum- 


bium 


but such a choice is facilitated by Table II, in 
which the metal to be evaporated and the 
filament materials are listed in the order of 


suitability. 
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bium 


In conclusion, the author acknowledges with 
pleasure the valuable suggestions which Dr. P. 
H. Carr offered during the progress of this 


investigation. 
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Steady-State Analysis of Multiconductor Transmission Lines 


Louis A. PiPEs 
Graduate School of Engineering, Harvard University. Cambridge, Massachusetts 


(Received May 1, 1941) 


Starting with the fundamental field equations of Maxwell, the various approximations 
involved in the usual engineering formulation of the multiconductor problem are examined. 
Subject to these approximations, the set of differential equations which is taken as the 
starting point of the classical transmission theory in the steady state are derived and the 
character of the various parameters involved in these equations is studied. By the use of the 
Laplacian transformation and certain theorems of matrix algebra these equations are solved 
for the case of general terminal conditions. Certain special cases where the system exhibits 


symmetry are discussed. 


INTRODUCTION 


HE mathematical theory of guided waves 

propagated along a system of uniform 
parallel conductors goes back to the early work 
of Heaviside and Kelvin. These investigators 
based their analysis on the concept that a trans- 
mission system can be represented as consisting 
of uniformly distributed coefficients of induct- 
ance, Capacitance, resistance, and leakage. This 
classical theory is based on the simple concepts of 
electric circuit theory and assumes that the 
circuit parameters are calculable from the geome- 
try of the system. 

However, by attacking the problem ab initio 
by the aid of the Maxwell field equations it is 
found that the form of the differential equations 
of the classical transmission theory is rigorously 
valid only for the ideal case of a system con- 
sisting of perfect conductors embedded in a 
perfect dielectric. This assumes no dissipation of 
energy in the conductors or in the dielectric. 
Nevertheless, since in any practically efficient 
system the amount of dissipation of energy is 
small, the distribution of fields and currents may 
be regarded as perturbed but slightly from the 
idéal dissipationless case. 

In this discussion the above program will be 
followed and it will be demonstrated that the 
electric and magnetic fields in the dielectric 
region and in the conducting region may be 
derived from two wave functions whose form 
may be deduced from the geometry and the 
electric constants of the system. From the nature 
of these wave functions, it is then shown that the 
complex currents and potentials of the system 
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satisfy the so-called ‘‘generalized telegraphist’s 
equations” which are the starting point of the 
classical or engineering treatment of the problem. 
By following this procedure, the physical signifi- 
cance of the circuit parameters of the classical 
formulation and their relation to the field equa- 
tions can be fully appreciated. 

After the form and general validity of the 
“generalized telegraphist’s equations’ has been 
established then by the use of matrix algebra and 
the Laplacian transformation these equations are 
solved and the complex current and potential 
distribution is determined to satisfy general 
boundary or terminal conditions. 


I. STATEMENT OF THE PROBLEM 


Let us consider a system of parallel, circular, 
cylindrical conductors embedded in a dielectric. 
The system may represent a group of overhead 
wires parallel to the surface of the earth, or a 
multiconductor cable enclosed in a sheath. The 
formal analysis will apply to both cases. However, 
it will be apparent that the calculation of the 
circuit parameters will be of greater difficulty for 
the cable because of the close juxtaposition of the 
conductors in this case. 

The n parallel conductors will be taken as 
parallel to the z axis of a Cartesian reference 
frame. It will be assumed that certain electro- 
motive forces are applied to the terminals of the 
conducting system at z=0, and that at z=S, the 
cylindrical conductors are connected together 
and to the ground or cable sheath by means of 
linear networks of known impedances. The as- 
sumption will be made that the length of the 
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cylindrical system is so great that the region of 
space in which the fields are perturbed by the 
terminals is negligible in comparison to the region 
between the terminals where this perturbation 
has no effect. 
The terminal electromotive forces will be 
assumed to be of a periodic nature of frequency 
f=w/2zx but of arbitrary phase. These terminal 
electromotive forces give rise to a _ periodic 
electric field E and a periodic magnetic field H in 
the space inside and outside the conducting 
system and also to a set of periodic currents 7, in 
the various cylindrical conductors. The problem 
under consideration is the computation of the 
field and current distribution from the prescribed 
electrical and geometrical data of the system. 


II. THe GENERAL FIELD EQUATIONS AND THE 
POTENTIALS 


The discussion will begin with the Maxwell 
field equations. This set of equations in a 
rationalized m.k.s. system of units takes the form 


VX E=—0B/dt, (1) 
VxH=a+0D/ot, (2) 
vV-B=0, (3) 
V-D=p, (4) 
B= ,H, (5) 
D=kE, (6) 
i=cE, (7) 


where, in this system of units, 


E =the electric intensity (volts/meter), 

H =the magnetic intensity (ampere turns/meter), 
ti =the current density (amperes/meter?), 

B =the magnetic flux density (webers/meter?), 

D =the electric displacement (coulombs/meter?), 
o =the conductivity (1/ohm meter), 

k, 1=the inductive capacities of the medium, 

k, =k/ko=dielectric constant or specific inductive capacity, 
Mr =u/uo=the permeability of the medium, 

ko= (1/367) X 10~° farad/meter, 

uo = 4X 1077 henry/meter. 


The factors wo and ko satisfy the fundamental 
relation: 


c=1/ (oko)? = 2.998 X 108 meters/second 


(velocity of lightin vacuum). (8) 
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The potentials 


The analysis of the electromagnetic field is 
facilitated by the use of auxiliary functions 
known as the potentials. This may be done by 
placing 

B=VxA, (9) 
E=—Vo—0A/adt, (10) 


where ¢=the scalar potential (volts), 


A=the vector potential (joules/ampere). 


In a medium of inductive capacities k and uy, 
and conductivity o, the potentials satisfy the 
wave equations 


V2 — kud?h/ dt? —ucdd/dt=0, 
V?A —kydA /dt?—ucdA /dt=0. 


(11) 
(12) 


The potentials are not independent but satisfy 
the relation: 


V-A+ypukd¢/dt+uco=0. (13) 
III. ELEMENTARY HARMONIC SOLUTIONS OF THE 
WAVE EQUATIONS 

For the purpose of discussing the propagation 
of periodic waves of frequency f=w/2m along a 
set of conductors parallel to the z axis it is 
convenient to assume solutions of the wave 
equations (11) and (12) of the form 


A=Ao(x,y)e™tivt, 


b= do(x,ye tim! 


(14) 
(15) 


where j= / —1, and a is a complex number to be 
determined. The complex number, a, is called the 
complex propagation constant. If these assumed 
forms of the solutions are substituted into the 
Eqs. (11) and (12), the following equations are 
obtained: 


V?Ao(x,y) = (jwou—a?—w?/v*)Ao, (16) 
V*g0(x,¥) = (jwou—a?—w*/v*)oo, — (17) 
where v is defined by: 
v=1/(uk)'=the velocity of light 
in the medium. (18) 
If we let 
b? = (jwon—w?/v’), (19) 
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then these equations may be written in the 
symmetric form 


V?Ao(x,y) = (b? —a?) Ao(x,y), (20) 


V*ho(x,y) = (b? —a*) do(x,y). (21) 


The problem of determining the periodic field 
and current distribution in a system of parallel 
conductors is reduced to the solution of (20) and 
(21) subject to the proper continuity and bound- 
ary conditions. It must be realized that the 
functions A o(x,y) and ¢o(x,y) are not independent 
but are related to each other through the 
keq. (13). 


IV. IpEAL CASE, PERFECT CONDUCTORS IN 
PERFECT DIELECTRICS 


Before discussing the general dissipative case, 
let us consider the ideal dissipationless case in 
which a set of parallel, perfect cylindrical con- 
ductors whose elements are parallel to the z axis 
are embedded in perfect dielectrics. In this case, 
because the dielectrics are assumed to be perfect 
insulators there cannot be any x or y components 
of current. It follows, therefore, that the cur- 
rents in the various conductors have only z 
components and hence the vector potential 
A,(x,y) has only a z component. That is, we have 


Ao(x,y) = (0, 0, A). (22) 


Since this vector has only one component there 
cannot be any confusion caused by denoting it by 
A (x,y) without the bar. We thus see that in this 
ideal case the fields are specified by two scalar 
functions A» and @». If we substitute (14) and 
(15) into (13) and make use of the relation (22), 
we obtain 


Ay=b' do jwa (23) 


as the relation between the two scalar functions 
Ay and ¢». To determine the fields, let us intro- 
duce the notation 


B= By(x,y)e-@tivt, (24) 
E= E,(x,y)e-=*tivt, (25) 

then from (9) and (10), we have 
By =idAo/dy —j0Ao/dx, (26) 
Ey = —iddy 0x —jdbo/dv+K(ado—jwAo), (27) 
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where i, j, and k, are unit vectors parallel to the 
x, y, and z axes. From this we see that there is no 
component of magnetic field parallel to the z axis. 
We now consider the form of the equations in the 
various regions. 


1. Dielectric region 

In the perfect dielectric, we have: 
(28) 
b? = —w?/y?. (29) 


Hence in the dielectric region, the functions ¢o 
and A,» satisfy: 


V*Ayp=- (a?+w? v*)Ao, 
V°*oo= — (a*+w? v*)do. 


(30) 
(31) 


Now in the dissipationless case under considera- 
tion, the energy flow must be in the z direction 
since no energy is being consumed in the perfect 
conductors or in the perfect dielectric, it thus 
follows that the Poynting’s vector, N=EXH 
(watts/meter?) must have only a z component. 
Because of this it is necessary that the electric 
intensity Ey cannot have a z component. From 
(27) this leads to 


(ado —jwA o) = (a—b?/a)oo=0 (32) 
in the dielectric. Hence, 
a? =)? = —y?/y?, (33) 


Therefore the propagation constant has the value 


a=+j(w/v). (34) 


From this we see that the propagation constant 
is an imaginary number in the dissipationless 
case. If f is the frequency of the various periodic 
quantities in cycles/sec. of the oscillatory electro- 
motive forces impressed at the terminals of the 
system, we see that approximately 


a=2f/3X108=2X10-*f (1/meter). (35) 


Hence in the ideal case, the propagation constant, 
a, is a pure imaginary and its absolute value is a 
very small number at all frequencies. As a 
consequence of (33), it is seen that (30) and (31) 
reduce to 


V*A,=0, (36) 


V°*oo=0. (37) 
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In view of (34), the relation (23) becomes 
Ao=¢0/V. (38) 


Hence in this ideal case, the Eqs. (26) and (27) 
reduce to 


Pag 
By =—-(iddo ‘OV —jddo Ox), (39) 
E,= —iddy Ox —jddo, ‘OY. (40) 
From this we obtain the relation 
E,- By, =0. (41) 


This shows that in the case of perfect con- 
ductors embedded in perfect dielectrics, the 
electric and magnetic vectors are both in planes 
normal to the z axis and normal to each other and 
in time phase. The flow of energy is parallel to the 
z axis. Since the propagation constant is given by 
(34) and is imaginary, we see that the periodic 
fields are propagated without attenuation and 
with the phase velocity v. 
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2. Conducting region 


In the conducting region, s—«, and the Eqs. 


(20) and (21) reduce to 
V?Ao(x,y) =jwopA, (42) 
V*oo0(x,y) =jJwoudo. (43) 
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These equations are of the form of the well- 
known heat diffusion equation. Since in this case 
we are assuming infinite conductivity, the solu- 
tion of these equations shows that there is no 
penetration of field into the conductors. Ac- 
cordingly, we may take ¢» and A» as constant 
over the conductor cross section. 

It is thus seen that in this ideal case there are 
no fields inside the perfect conductors and that 
the fields outside the conductors are given by 
(39) and (40). Both the electric and magnetic 
fields are derivable from a single scalar potential 
function, ¢9 that satisfies Laplace’s Eq. (37) in 
the dielectric region and reduces to prescribed 
constant values at the conductor cross sections. 
From these considerations and the form of the 
propagation constant as given by (34) it is seen 
that in the ideal case the conductors guide the 
electromagnetic waves without distortion or 
attenuation. The waves travel with the phase 
velocity : 


v=1/(uk)'=c/(u,k,), (44) 


where y, and k, are the relative permeability and 
dielectric constant of the insulating medium and 
c is given by (9). 


3. The engineering formulation of the ideal case 


Let us now consider how the above field 
analysis leads to the formulation of the problem 
in the usual classical circuit or engineering 
manner. To fix the ideas, let us consider the 
problem of parallel cylindrical conductors over 
a parallel earth plane as shown in Fig. 1. In this 
ideal case, we consider the conductivity of the 
conductors and of the earth as being infinite. 
From the above considerations, it is seen that the 
electric and magnetic fields may be specified by 
the single function ¢0(x,y). This function satisfies 
Laplace’s Eq. (37) and assumes constant values 
over the cross section of the conductors and the 
earth. It is thus evident that the determination 
of the field distribution is determined from the 
potential or electrostatic problem of obtaining a 
harmonic function @¢9 that reduces to various 
prescribed constant values throughout the con- 
ductor cross sections and at the earth’s surface. 

(a) Instantaneous and complex quantities.—To 
facilitate the discussion and to arrive at the 
engineering formulation of the problem we intro- 
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duce the following notation. Let the instantane- 
ous charge per unit length on conductor r be 


q,(z,t) = ore tivt = O,eiv', (45) 


Let the instantaneous potential at the point z of 
the rth conductor be 


$,(z,t) = ho,e-e tivt = Veit, (46) 


Let the instantaneous current flowing in the rth 
conductor at the point z be 


i,(z,t) _ Te +jwt — Tei (47) 


then the quantities Q,, V, and J, are the complex 
charge per unit length, the complex potential and 
the complex current used in conventional alter- 
nating current circuit theory. 

Let us take the potential of the perfectly con- 
ducting ground plane as being zero. At any point 
z, let the instantaneous potentials relative to the 
ground plane be ¢,(z,t) throughout the cross 
section of the rth conductor. The determination 
of the electric and magnetic fields from the 
potential distribution depends on the solution of 
the well-known electrostatic problem where it is 
required todetermine a harmonic function ¢o(x,y) 
which reduces to zero at all points of the ideal 
ground plane and to the constant value do, 
throughout the cross section of the rth conductor. 
The solution of this problem may be most easily 
accomplished in the case of Fig. 1 by the method 
of images. In any case, the potential of the rth 
conductor may be expressed in the form 


Por = Prigort proqo2t *** +Prngon, (48) 


where the p,, quantities are the Maxwell coeffi- 
cients of potential. 

(b) The equation of conservation of electricity.— 
In order to obtain the ordinary classical or 
engineering formulation of the multiconductor 
problem, it is necessary to write the equation of 
conservation of electricity for the case under 
consideration. If we take the divergence of both 
sides of Eq. (2), we obtain 


V:-(VXH)=0=V-(a+0D/at). (49) 


In a medium of conductivity o and inductive 
capacity, k, we have: 


V-(o+kd/at) E=0. (50) 
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If we now apply Gauss’ theorem to the surface 
bounding a length dz of the rth cylindrical con- 
ductor as shown in Fig. 2, we obtain: 


ate, = = $ 0x36. ‘On)ds, (51) 


where io, is the amplitude of the current flowing in 








Fic. 2. 


the metal. The directions m and s are normal and 
tangential to the conductor as shown in Fig. 2. 
The quantity ¢,, is the complex conductivity of 
the dielectric medium surrounding the rth con- 
ductor. It is defined by: 


o-=oatjwka, (52) 


where o, is the conductivity and k, is the in- 
ductive capacity of the dielectric. In the ideal 
case under consideration, ¢g=0 and (51) reduces 
to: 


aio, = —jwka ¢ (Ago,/dn)ds. (53) 


By an application of Gauss’ theorem to Eq. (4) 
over the same surface of Fig. 2, we have: 


dor = -bag (860, ‘On)ds. (54) 


Hence (53) reduces to: 
ior =jWQor- (55) 


Now the set of m linear equations, (48), may be 
solved for the charge amplitudes in terms of the 
potential amplitudes and written in the form 


Jor =CripoitCropo2t * * * +CrnGon, (56) 


where the c,, quantities are the coefficients of 
capacity and induction. Therefore (55) may be 
written in the form: 


Ato, =jw(CridoitCropo02+ ill +Crnon)- (57) 
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Then by (46) and (47), we obtain 
—dI,, ‘dz =jw(Cr Vitcre Vot ~ +CrmnVn), 


r=1,2,---n. 


(58) 


The set of Eqs. (58) is the conventional set 
which determines the variations of the complex 
conductor currents with distance in the ideal non- 
dissipative case. 

To determine the variation of potential with 
distance, let us begin with the Eq. (32) expressing 
the vanishing of the z component of Eo. We 
therefore write 


doy =jwAo. (59) 
In view of (38), we therefore have at the 
surface of the rth conductor: 


Ador = jwoor/v. (60) 
By substituting (48), we may write (60) in the 
form 


ado, =—(Pridoit Pr2go2+ iit + PrnGon)- 


Vv 


(61) 


However, from 


relation: 


(55) and (34) we have the 


Jor = Q10,/jw = 10/9. (62) 
Hence we may write (61) in the form 
.... ; 
Abo, =—(britort prato2 t+ i is + Prnton). (63) 


In view of (46) and (47) this gives 


jw 
~(Prlit Prolet ss + Prat n). (64) 





—dV,/dz= 


¢ 


This set of equations gives the variation of the 
complex potentials and currents of the system. 
It is customary in the classical literature to write 
this set of equations in the form 


—dV,, ‘dz =jw(LilitLrele+ Pee +L,.I,), (65) 


r=1, > ++ 9M, 


where the coefficients L,, are called the coeffi- 
cients of self and mutual inductance of the lines. 
Comparing (65) and (66), we thus see that the 
coefficients of inductance are related to the 
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coefficients of potential by the equation 


Lis= prs/v? henrys/meter. (66) 


The relation (66) holds in general for the ideal 
dissipationless case for all conductor sizes and 
spacing. The coefficients of potential, p,,, may be 
calculated from the geometrical arrangement of 
the system by the usual methods of potential 
theory. In the ideal cable, the sheath replaces the 
ground as the perfect conductor at zero potential 
and the above analysis is applicable. 

From the above discussion, it is evident that 
the classical or engineering treatment is com- 
pletely rigorous for the dissipationless case. In 
such a case there is only one velocity of propa- 
gation, that of light in the dielectric medium. 
The system may then be completely specified in 
terms of its coefficients of potential. 


V. THe CASE OF FINITE CONDUCTIVITY 


We now turn to the consideration of the actual 
case encountered in practice. The conductors will 
no longer be considered perfect nor will the 
dielectric be considered as perfectly insulating. 
This actual case will be considered as a sort of 
perturbed case of the ideal one. It will be neces- 
sary to make certain approximations which will 
now be enumerated. 


1. General approximations 


(a) It will be assumed that the propagation 
constant, a, is an extremely small quantity whose 
absolute value is given approximately by (34). It 
will be assumed that the real part of a is small 
compared to its imaginary part. Since |a| is of 
the order of magnitude of 10-* f/meter where f is 
the frequency of the periodic electromotive 
forces, it is seen that a is a very small quantity 
even for frequencies of several millions of cycles 
per second. 

(b) The conductivity, o, of the conductors 
will be assumed to be of the same order of 
magnitude as that of copper (5 X 107 mho/m) and 
the relative dielectric constant k, of all materials 
will be taken as being of the order of unity. 

(c) The conductivity of the insulators will 
be taken to be of the order of (10-* to 10~'° 
mho/meter). 

(d) The conductivities of geological materials 
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vary greatly from place to place. However, the 
conductivity of wet ground is of the order of 
magnitude of 10-? to 10-* mho/meter and of dry 
ground between 10~-* to 10-5 mho/meter. Since 
the conductivity of most dielectrics is of the order 
of 10-* to 10~-'® mho/meter, it is evident that the 
ground may be considered as a good conductor. 

(e) The relative permeability 4, will be as- 
sumed to have the value 1 in all dielectrics and 
non-magnetic materials. 

(f) Because of the low conductivity of the 
dielectric medium, the component of the mag- 
netic field /7, in the z direction will be everywhere 
infinitesimal since it could only be produced by 
the x—y components of current which are 
infinitesimal in the dielectric. Therefore, as in the 
ideal case of perfect conductors and insulators the 
magnetic vector potential, A, will be assumed to 
have only the z component, A,. 


2. The general equations 


Because of the above assumptions, it is seen 
that the problem may be regarded as a pertur- 
bation of the ideal case discussed in Part I. 
Accordingly, we return to the Eqs. (20) and (21) 
and consider the form that these equations take 
in the dielectric and conducting regions. 

(a) Dielectric region.—In the dielectric region, 
we have 


V?Ao(x,y) = (bu? —a?)Ao(x,y), (67) 
V2h0(x,¥) = (ba? —a*)ho(x,¥), (68) 


where 
ba? = (jwoapa— w/v") (69) 


and o, and yw, are the conductivity and perme- 
ability of the dielectric region. Now in the 
dielectric, b,* and a? are very small quantities and 
are nearly equal. We may, therefore, to a high 
degree of accuracy, place (b,?—a*)=0. 

To this degree of approximation, therefore, 
the functions ¢@o9 and Ay» satisfy the Laplace 
equations 
, V°6o=0, (70) 


V*A,o=0, (71) 


in the dielectric region. If we now apply the 
condition V-E=0, to the electric intensity in the 
dielectric region, we obtain 


— aE o,= + (0*b0/dx?+07%o0/dy?) =0. (72) 
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This shows that in the dielectric region this 
above approximation leads to the condition that: 


(electric transverse field in dielectric region) 
> (longitudinal field). 


This is consistent with the consideration that 
the energy flow as expressed by the Poynting’s 
vector 


N=E XA (watts/meter?) (73) 


shall have its longitudinal component much 
larger than its transverse component. 

(b) The conducting region.—As a consequence 
of the assumption that H,=0, it follows from the 
Maxwell equation (1) that in the conductors the 
electric transverse field is of the order of magni- 
tude of a/woy compared to the longitudinal 
electric field, Eo.. Since this is a very small 
quantity, we are led to the condition 


(electric transverse field in the conductor region) 
«(longitudinal field, E.0). 


From an energy consideration, this is consistent 
with the assumption that the dissipation of 
energy in the conductors due to longitudinal 
currents is very large compared to that of 
charging currents. From this, it follows that 


Voo=0 (74) 


throughout the conductor cross sections, and 
hence the function ¢@o9 assumes the constant value 
go, throughout the cross section of the rth 
conductor. Since a is a very small quantity, we 
may take 


Fo: = —JwAy (75) 


inside the conductors. In view of (26) this gives 


. 1 
By=—(—idEo, dy+joEo. Ox) (76) 
jw 


for the magnetic induction inside the conductors. 
Since b?=jwop—w*/v? and a* and w*/v? are small 
quantities of the same order of magnitude, the 
magnetic vector potential, Ao, satisfies the 
equation: 


V?A o(x,y) =jwopA o. (77) 


This is the well-known heat diffusion or ‘‘skin- 
effect’’ equation. Because of the relation (75), the 
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function, Eo. also satisfies the equation 


V°Eo, =jwopEy,. (78) 


The quantities o and yp in (77) and (78) are the 
conductivity and magnetic inductive capacity of 
the conducting medium. 


3. The continuity or boundary conditions 


The analytical solution of the multiconductor 
problem depends upon obtaining solutions of the 
various equations in the different regions and 
then adjusting these solutions to satisfy the 
various continuity or boundary conditions. These 
will now be listed. 

(a) The longitudinal electric intensity, Eoz.— 
The longitudinal component of electric intensity, 
E», is of finite order of magnitude everywhere and 
is continuous at all boundaries. Since E»o,=ado 
—jwAy, it follows that A» and aq» are also finite. 
If @) and A» are continuous at all boundaries, 
then Eo. will be continuous. 

(b) The magnetic vectors—Let us consider the 
cross section of the rth conductor as shown in 
Fig. 2. Then if ” is the normal and s the tangential 
direction along the conductor cross section, we 
have, as a consequence of (9), the relations 


Bos; = —0Ao/On, 
Byn=0A 0 ‘OS. 


(79) 
(80) 


Now at the boundary between two materials of 
different magnetic permeability, it is a conse- 
quence of the field equations that the normal 
component of B and the tangential component of 
H be continuous. As a consequence of (80), it is 
seen that the continuity of Bo, is assured by the 
continuity of Ao. To assure the continuity of Ho,, 
it is necessary from (79) that 


ii ‘On 
bu 
be continuous at all boundaries. 

(c) The continuity of the total current.—By an 
application of Gauss’ theorem to the Eq. (49) 
expressing the continuity of total current to the 
boundary between two different dielectric re- 
gions, it is evident that the quantity 


g.0go/dn, 


where o, is the complex conductivity defined by 
(52) must be continuous at these boundaries. 
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(d) Relations involving the conductor currents.— 
As a consequence of (2) and the large value of the 
conductivity assumed for the conductors we have 
the relation 

—1 
lor oo (0A Or ‘On)ds, 
be 


(81) 


where the integration is performed in the 
tangential, s direction at the surface of the rth 
conductor. If we integrate the ohmic relation (7) 


across the normal area S, of the rth conductor, we 


obtain 
ino { EadS, 


where o, is the conductivity of the rth conductor. 


(82) 


4. Analytical formulation of the problem 

From the above general analysis, it is seen that 
the dissipative multiconductor problem may be 
formulated in terms of three scalar functions, 
oo, Ao and Eo, that must satisfy the following 
conditions: 

(A). The function oo 


(a) V2¢o=0 everywhere. 

(b) ¢9=¢0, throughout the cross section of the 
rth conductor. 

(c) $0 is continuous everywhere. 

(d) o.0¢0/0n is continuous at the boundaries be- 
tween different dielectrics. 

(e) @o=0 at the surface of the ground plane or 
the conductor sheath. 

(f) go-=—kRG (Ado,/dn)ds at the surface of the 
rth conductor. 


(B). The functions Ay and Eo, 
(a) V*Ao=0 in dielectrics. 


(b) Ao is continuous at all boundaries. 


1 
(c) —dAy/dn is continuous at all boundaries. 
m 
. ee 1 ; 
(d) io-=—— ® (0Ao,/dn)ds at the surface of the 
m 


rth conductor. 


(e) V?Eo.=jwo,u-Eo, inside the rth conductor. 


(f) in=or| f EudS, across the rth conductor. 
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Since the magnetic induction vector may be 
expressed by (76) in the region inside the con- 
ductors and by (26) in the region outside the 
conductor, the continuity conditions of H, and 
B, at the surface of the rth conductor may be 
used to obtain the relations 


—- 

(g) (0A 9 on), = : (OE, On),, 
jwu, 

(h) (0A 0 OSs), = —(dE». OS)». 
jw 


The relation (g) expresses the continuity of the 
tangential component of 7 at the surface of the 
rth conductor, uw, is the permeability of the 
conductor and yu that of the dielectric. Equation 
(h) expresses the continuity of the normal com- 
ponent of Bo at the surface of the rth conductor. 


5. The engineering formulation of the dissipa- 
tive case 


(a) The coefficients of capacity and leakage.— 
We now turn to the engineering or classical 
formulation of the dissipative multiconductor 
problem. It will be found that the solution of the 
above equations for 0, Ao, and Eo, may be 
formulated in terms of the well-known parame- 
ters of the classical or engineering theory. Let us 
first consider the differential equation and the 
boundary conditions satisfied by the scalar po- 
tential function ¢o(x,y) listed in Section (A) 
above. We see from these equations that the 
determination of this function when the geome- 
try of the conductors and insulators is specified is 
the solution of the well-known two-dimensional 
potential problem. We may, therefore, as in the 
ideal case write 


dor = Prido1t+ Pr2do2+ oe + PrnQons (83) 


where ¢o, is the potential of the conductor r. The 
brs Quantities are the Maxwell coefficients of 
potential and the go, quantities are the various 
charge amplitudes on the conductors. Let us now 
consider the Eq. (51) expressing the continuity of 
the total current along the rth conductor. This 
equation may be written in the form 


Cig = — g (cat+jwks)(Ad0,/dn)ds. (84) 





The quantities ¢, and k, are the conductivity and 
the specific inductive capacity of the insulating 
medium surrounding the rth conductor. Now 


tr = § odo ‘On)ds (85) 


is the leakage current per unit length from the 
surface of the rth conductor. By the use of Eq. 
(54) giving the relation between the charge per 
unit length and the normal derivative of the 
potential we may write (84) in the form 


Ao, = (jwteoa, Ra) Qor- (86) 


Since as in (56) we may express the charges in 
terms of the coefficients of capacity and induction 
and the potentials, we may write 


a lor — (jwt+oa/ka) 
x (Crido1 +Cr2602+ ‘shes +Crnon)- (87) 


In the analysis of open wire lines, it is custom- 
ary to write the Eq. (87) in terms of leakage 
coefficients. This may be done by using the 
notation 


£,r=leakage conductance to ground in mhos/ 
meter, 

£,;=leakage conductance between conductors r 
and s in mhos/m. 


In terms of these coefficients, the leakage con- 
ductance to ground and to the other conductors 
from’ the rth conductor may be expressed as 
follows 


1 r= £r1( bor — Go1) 
+ 2,2(b0,— $02) + elie + 2rror 
= hor(Zrtget inh + 2rn) 
— 271901 — 2,202 7? *ZrnPon- (88) 
If we now let 
Grr =2rntZet pei +2rn; 
(89) 
G,.=G,,= —2rs = — Ler, 


we may write the leakage current of the rth con- 
ductor in the symmetrical form 


trr=Grib01t+Grodo2+ +» +Grndon: (90) 


From this it is seen that we may write the Eq. 
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(87) in the form 
ator = (GridoitG,2po2+ * > - +Grndon) 
+jw(Cridoitcrapo2+ + +++Crnbon)- (91) 
By (46) and (47), we may write 
—dI,/dz= VinVit Vr2Vet-- ++ ¥inVn, 
(r=1,2,---m), (92) 


where 


Vs=Gret+jwers. (93) 


The set of Eqs. (92) is the conventional set 
giving the variation of the complex currents, J, in 
terms of the potentials and admittances of the 
system. A form of the equations which is more 
applicable to the case of cables where the con- 
ductors are embedded in homogeneous insulators 
may be obtained from (87) is the following one 


—dI,/dz=(jw+oa/ka) 
X (Cn VitceVet+: , >+OrnV2). (94) 


It will be noted that Eq. (92) is similar in form 
to Eq. (58) for the ideal case. The only difference 
is the addition of the leakage coefficient terms in 
the dissipative case. 

(b) The coefficients of resistance and inductance. 

-We now turn to a consideration of the form 
that the Eqs. (65) take in the dissipative case 
under consideration. From the above analysis, we 
have 

Eo: = ado —jwAy (95) 


in the region outside conductors, while Eo, 
satisfies the conditions (e), (f), (g) and (h) of 
Section 4, (B) inside and on the boundary of the 
conductors. There are two cases to be considered. 

1. The case of negligible proximity effect: If 
we assume that the separation of the conductors 
from each other and from the earth is so great 
that the disturbing effect on the field of the rth 
conductor produced by the fields of the currents 
in the other conductors and the earth may be 
taken as negligible, then we have negligible 
‘proximity effect.’’ This is the usual case in open 
wire lines. In this case, the Eq. (e) satisfied by Eo. 
in the rth conductor is 


1 
d?E,/dp?+—dEo./dp=jwoou,Eo:, (96) 
p 
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where p is the radius vector of Fig. 3 drawn from 
the center of the rth conductor. Since proximity 
effect is being neglected, it is assumed that Eo, is 
independent of 6. If the rth conductor is a solid 
wire of radius a,, the integration of (96) subject 
to (f) gives 








Ever =1orZrr (97) 
A 
qa |F 
9 — 
Fic. 3. 


where 2,, is the internal impedance of wire r and 
is of the form 


Serr = Pert JW rr. (98) 
If we denote by 
¥o9=1/0,ra,? (ohm/meter) (99) 
the d.c. resistance of the rth conductor, then 


Zrr =rona,J o(na,)/2J,(na,), (100) 
where ae 
n=jr/j(o,Wpr)! (101) 


and Jo is the Bessel function of the first kind and 
zeroth order while J; is the Bessel function of the 
first kind and first order. The real part of z,, is the 
so-called skin-effect resistance of the rth con- 
ductor, 7,, and the imaginary part is the internal 
inductance of the wire, I,,. 

The equations which the function Ao satisfies 
serve to determine the quantity Apo as a linear 
function of the m conductor currents in the form 


AveE flea en (102) 


The quantities f,(x,y) are functions depending 
on the geometry of the conducting system and on 
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the electrical constants. Substituting (97) and 
(102) into (95) when evaluated at the surface of 
the rth conductor, we obtain an expression of the 
form: 


n 
Srrtor = Ado; — : Zrktor- (103) 
k=1 

This is a generalization of the Eq. (61) of the 
ideal case. By comparison of these two equations, 
it is seen that in the ideal case the internal 
impedance 2z,, is zero and the impedance functions 
Z,, reduce to pure imaginary terms of the form 
jwL, where L, are the coefficients of mutual 
inductance between the conductors when r#k 
and the coefficient of external self-inductance of 
conductor r if r=k. In the dissipative case under 
consideration, however, the impedance functions 

Z,, are complex and of the form: 


Lo = RR. +jwL. (104) 


If the conductors are small and widely separated 
from each other and from the earth, then the 
resistance terms R,, become negligible. In such a 
case the coefficients L,, approximate the values 
given by (66) of the ideal case. By (46) and (47) 
we may write (103) in the form: 


n 

—dV,/dz=2,-1,+>, Zul (105) 
k=1 

in terms of the ordinary complex potentials and 

currents of the system. 

2. The case of appreciable proximity effect: 
In cable circuits, at telephonic frequencies be- 
cause of the close juxtaposition of the conductors, 
the disturbing effect of the field in the rth by the 
fields of the currents in the neighboring con- 
ductors may not be neglected. In that case, the 
electric field Eo, in the rth conductor does not 
satisfy Eq. (96) because Eo, is no longer inde- 
pendent of 6. In this case the equation satisfied 
by Eo, is 


1 0 OF. 1 07Eo, ; 
~ —(» “) += ——=Jwo,p,E:. (106) 
p Op Op p* 00 


The solution of this equation subject to the 
boundary conditions of Section 4, (B) may be 
written in the form 


Eon = > Zrklok, (107) 


k=1 
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where Ep,, is Eo, evaluated at the surface of the 
rth conductor. The additional terms arise from 
the proximity effect. In this case we may write 
(105) in the form 


n 


—dV, dz=>> 2nd +> belli (108) 
k=1 k=1 
It is customary to incorporate the self and mutual 
internal impedance coefficients 2z,, with the 
external impedance coefficients, Z,, and write 
(108) in the form 


n 
—dV, dzs=>° “rid k, (109) 
k=1 
where in (109), the coefficients Z,, include the 
internal impedance and proximity effect coeff- 
cients as well as the external mutual impedance 
coefficients. 

We have thus arrived at the form of the 
equations giving the variation of complex po- 
tentials and currents in the dissipative multicon- 
ductor lines. The actual calculation of the 
impedance parameters from the geometry and 
electrical properties of the system is usually a 
matter of great difficulty. These parameters are, 
in any case, experimentally determinable. 


VI. THE GENERALIZED TELEGRAPH EQUATIONS 
AND THEIR SOLUTION 


1. General considerations 


It is seen from the above that the analysis of 
multiconductor transmission lines may be formu- 
lated in terms of the two sets of Eqs. (92) and 
(109). These equations are the ones used in the 
elementary or classical transmission theory and 
we see that they are consistent with the field 
equations to the approximations involved, that 
is, that the actual system be slightly perturbed 
from its ideal dissipationless state. We may 
write these equations in the matrix form 


d 

—-—(U)=LY](V), (110) 
dz 
d 

—-—(V)=[Z])), (111) 
dz 


where (J) and (V) are column matrices whose n 
elements are the complex currents and potentials 
to ground or the cable sheath of the conductors, 
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respectively. [Y] and [Z] are non-singular 
square matrices of the mth order whose elements 
are the admittance and impedance coefficients 
Y,, and Z,, defined by (93) and (109), respect- 
ively. These equations are called the ‘‘generalized 
telegraph equations” in the literature. By a 
formulation of the energy distribution and dissi- 
pation of the system it may be demonstrated that 
the elements of the matrices [ Y] and [Z] satisfy 
the relations: 


Vis= Yor, (112) 
LZrs=ZLers (113) 


and hence the square matrices, [ Y ] and [Z ] are 
symmetric. We have thus seen that the Maxwell 
field theory leads to the engineering or classical 
form of the equations which govern the propa- 
gation of the complex currents and potentials of 
the multiconductor lines in the region that is so 
far removed from the terminals of the lines that 
the perturbing effect of the terminal apparatus 
may be neglected. 

We now make the assumption that the Eqs. 
(110) and (111) continue to hold throughout the 
entire length of the conductors including the 
terminals of the system. It will be assumed that 
at one terminal which will be designated the 
“generator end’”’ there are certain devices that 
maintain alternating potential differences be- 
tween the various conductors and the ground or 
cable sheath of the same frequency f=w/2z but 
of different magnitudes and phases. At the other 
terminal of the multiconductor lines, it will be 
assumed that the conductors are connected to 
each other and to the ground or sheath by means 
of a certain system of conductors whose self and 
mutual impedances are known. This terminal 
will be called the receiving end. The problem of 
determining the potential and current distri- 
bution in this system is then one of solving Eqs. 
(110) and (111) in terms of the terminal boundary 
conditions. 


2. Solution of the equations, boundary conditions 


Let us consider the schematic diagram of Fig. 4. 
The elements of the matrices (V), and (J), are 
the complex potentials to ground and currents at 
the generator end of the various conductors of the 
system. The elements of the matrices (V), and 
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(I), are the complex amplitudes of the potentials 
and currents at the receiving end. In Fig. 4, the 
arrangement has been drawn as if it were a 
single conductor line. This is done for con- 
venience and it must be realized that the various 
matrices have n elements. The elements of the 
square matrix [Z ], represent the various imped- 
ance elements of the several lines to each other 
and to ground. In general, we have the relation 


(V),=[Z](Dr. (114) 


We must solve the Eqs. (110) and (111) 
subject to the boundary conditions at the 
terminals. In certain problems of power trans- 
mission, the receiving end currents and potentials 




















iI, 
~TII] 4 #1, 
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| cuba nea 
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are specified and it is required to determine the 
generator end currents and potentials. In other 
problems, we are given the generator end po- 
tentials and currents and the receiving end 
currents and potentials must be specified. In 
still other problems the potentials at the gener- 
ator end are specified and the receiving end 
impedances are given. In this case the current 
and potential distribution must be determined. 

The Eqs. (110) and (111) are two sets of linear 
differential equations of the first order with 
constant coefficients. The simplest method of 
solution is to introduce a Laplacian transforma- 
tion with respect to the space variable, z. That is, 
we introduce two column matrices, (£),; and (J), 
which satisfy the relations 


(V).= f e->+(V) dz, (115) 


(D.=p J e-*(T)dz. (116) 
0 


That is, every element of the matrices (£), and 
(J), are the Laplacian transforms of the elements 
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of the matrices (Z) and (J). Following the nota- 
tion of (5), the relations (115) and (116) will be 
expressed symbolically in the form 


(V):=(V), (117) 
(1).= (J). (118) 


From the basic theorem involving the transforms 
of first derivatives, we write 


d 
: (V)=—p(V),+p(V)., (119) 
az 
d 
P (1)= —p(),+p()., (120) 


where (V), and (J), are matrices whose elements 
are the values of the variables at z=0. By the use 
of the relations, (110) and (111) transform to 


P(V),—P(V) :=(Z (Ds, (121) 
PUD) ,—p(D:=LY](V)«. (122) 


Eliminating (£), from these equations, we 
obtain 


(P?U-LYIZ)(D).=p*D.— PLY MV). (123) 
Eliminating, (J);, we have 
(pP?U—-[ZILV))(V):=p?(V),—PLZ]D)o, (124) 


where U is the unit matrix of the mth order. Let 


(kJ=LYILZ] (125) 

and 
p?*=m, (126) 
[f(m) 1=[mU—k]. (127) 


With this notation, we have (123) in the form 
Lf(m)\D).=p?(Do—-PLY\(V),. (128) 
We thus have explicitly 
(1).=Li(m) T'Le?(Do—PLYM(V),]. (129) 
In order to compute the matrix [f(m)]", 
we may use the identity (see reference 6, 
Bibliography) 
[f(m) }-'=(F(m) ]/D(m), (130) 


where [ F(m) ] is the adjoint matrix of [f(m) ] and 
D(m) is the determinant of the matrix [f(m) ]. 
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Now if we assume that the determinantal 
equation, 


D(m)=0, (131) 


has n distinct roots, (m1, m2, ---m,) and the 
degree of the adjoint matrix [ F(m) ] is less than 
that of the polynomial D(m), we have 


n [ F(m,) ] 1 
f(s) F'*= 2 ———_ ————_.._ (132) 
r=1 D’(m,) (m—m,) 
where 
d 
D’(m,) =——D(m) |m =m, (133) 


dm 


Substituting this result into (129), we obtain: 


(1).=>- [F(m,) ]p?(1),/D’(m,)(p?—m,) 
r=1 


a [F(m,) |LY ](E),p/D’(m,)(p?—m,). (134) 


To obtain the complex amplitude matrix, (J), 
we must compute the inverse transform of (J),. 
This may easily be done by using the transforms 


p?/ (p?—m,)= cosh (/m,2), (135) 
p/(p?—m,)= sinh (\/m,z)//m,. (136) 


Substituting these values into (134), we have: 


(1) => [F(m,) }(J), cosh (4/m,z) /D’(m,) 


r=1 
— > LF(m,) LY \(E), 
r=1 
Xsinh (4/m,z)/D'(m,)\/m,. (137) 


This is a formal solution for the complex 
current amplitudes along the line in terms of the 
generator currents and potentials. If we let 
a?=m, d,-=++/m,, we see that in this case we 
have 2n propagation constants a,, which merge 
together to give the solution (137). If we have 
only a single line and its return, then the solution 
(137) reduces to: 


I=I, cosh (az)— V,(Y/Z)' sinh (az), (138) 


where a=\/ZY. 
To obtain the potential distribution, we must 
find the inverse transform of (V), of Eq. (124). 


Let 
(wJ=(Z)LY] (139) 
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and 
[g(m) ]=(mU—w]. 


It is easy to show that the determinant of the 
matrix [g(m)] is the same as the determinant 
D(m) of [f(m) ] above. This follows from the fact 
that the matrices [Y] and [Z] are symmetric 
matrices. We thus have the relation 


(140) 


| g(m) | =D(m). (141) 


Hence the roots of the determinantal equation 
g(m) | =0 are the same as those of the equation 
D(m)=0. We may now expand the matrix 
[g(m) | in the form 


[e(m) "=> [G(m,)/D'(m,)(m—m,), (142) 


where [G(m,)] is the adjoint of the matrix 
[g(m) |] evaluated at m=m,. Substituting the 
above relations into Eq. (124) and making use of 
the transforms (135) and (136), we obtain: 


(VY) = > [G(m,) ](V), cosh (24/m,) /D’(m,) 


-¥ [G(m,) IZ), 


Xsinh (24/m,)/D'(m,)/m,. (143) 


This is a formal solution for the complex 
potentials in terms of the generator currents and 
potentials. For the case of a single wire and its 
return circuit, this reduces to: 


V=V, cosh (az)—I,(Z/Y)' sinh (az). (144) 


3. The determinantal equation 


The Eqs. (137) and (143) give the formal 
solution for the current and potential distribution 
along the multiconductor lines in terms of the 
sending end conditions provided that the de- 
terminantal equation 

D(m) = |mU—-LY ][Z]| =0 (145) 
has n distinct roots, (m1, m2, ---m,). In this case 
we have the 2” propagation constants, ++/m,, 
y=1,2, ---n., 

In certain cases of practical importance, it will 
be found that the determinantal equation (145) 
has multiple roots and that the number of 
distinct propagation constants is reduced. We 
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will discuss the more common cases that occur in 
practice. 

(a) Ideal case, perfect conductors, perfect ground 
conductivity.—This case was examined from the 
field point of view in Section IV. Let us now 
consider it from the point of view of the ‘“‘gener- 
alized telegraphist’s’’ equations. Since for this 
case we have no energy loss, we have the 
relations 


[G]=[0], (146) 
[R]=[0], (147) 
LY]=julc], (148) 
[Z]=ju[L], (149) 
LY]LZ]=—-w*[cl[L}. (150) 


From Eqs. (48) and (56), we obtain the 


relation 
[eJ=[p}° 


where [p ] is a square matrix whose elements are 
the various coefficients of potential. For this ideal 
dissipationless case, we have 


[LJ=[p]/2° 


as given in Eq. (66). 
Substituting (151) and (152) into (150) one 
obtains 


CYJ(Z)=—w*[p tL J/v? = —w?U/2’, 


where U is the unit matrix of the nth order. If 
this value of [ Y ][Z] be substituted into (145), 


we have 


(151) 


(152) 


(153) 


D(m) = | (m+w?/v®) U| =0. (154) 
Or 


(m+w?/v?)"=0. (155) 


In this case the determinantal equation has n 
equal roots, 


= —wy?/y? (156) 


and there are only two propagation constants 
+jw/v, one associated with waves traveling in 
the +2 direction and one with waves traveling in 
the —z direction. In this case, it is easier to 
obtain the current distribution directly from the 
Eq. (123) rather than from the general equation 


(137). Substituting (153) into (123), one obtains 
(p?U+w*?U/v*)(1).=p?(Do— PLY (E)o. (157) 
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Or 
(1) 1:= (1) op?/ (p?+-w?/v?) 
— p[ Y](E),/(p?+w?/v?). (158) 


Using the transforms 
p*/ (p?+a*)=cos (az), (159) 
b/ (p?+a*)= sin (az) /a, (160) 
we obtain 
(1) = (J), cos (wz/v) 
—vl_ Y](E), sin (wz/v)/w. (161) 


In the same way, we obtain the potential 
distribution 


(V)=(V), cos (wz/v) 
—v[Z (J), sin (w2/v)/w. (162) 


(b) Transposed conductors, ring symmetry.—In 
many cases of practical importance, the con- 
ductors are identical and are arranged or trans- 
posed in such a manner that the matrices [Z ] and 
[Y] exhibit so-called ‘“ring-symmetry.”’ As a 
special case, let us consider the arrangement of 
the six-conductor cable of Fig. 5. This figure 
represents the cross section of a six-conductor 
cable whose conductors are arranged in a sym- 
metrical manner inside the cable sheath. If we 
assume that the conductors of the cable are 
identical and symmetrically spaced with respect 
to each other and to the cable sheath, then the 
impedance matrix of the system has the form 


'Z Zu Zs Zu Zs 211] 
Zi2 Z Zi2 Z13 Zi4 Z13 
r Zi13 Zi2 & Zi2 Zi3 Zi4 

Zizi, ‘ “ - , 163 
2i"\2, Zu Za £ Zu Zui °™ 
Zi3 Zi4 Zi3 Zi2 Z Zi 
| Z12 Zi3 Zis Zi3 Zi2 Z 








Here Z is the self-impedance of the various 
identical conductors and the mutual impedances 
between conductors at equal distances from each 
other around the ring are set equal to each other. 
Because of the symmetry in this case, the mutual 
impedance coefficients are only three in number. 
The admittance matrix of the system, [ Y ], has 
the same type of symmetry as [Z]. It may be 
demonstrated that if the matrices [Z] and [Y] 


796 


have the above symmetry then the determinantal 
equation (145) has m distinct roots. The im- 
portance of the ring-symmetric case is that it 
may be solved explicitly and is of practical 


S 
S § 
S _§ 

S 


Fic. 5. 


importance in the case of cables and transposed 
multiconductor aerial lines. Rather than to 
attempt to obtain the various roots of the 
determinantal equation directly, the Eqs. (121) 
and (122) may be transformed by means of the 
non-singular symmetric matrix, [A ], defined by 


[A ]=[4.], (164) 

where (r=1, 2, ---m), (s=1, 2, ---m) and 
A,=a- YD), (165) 
a=e", j=/-1. (166) 


This transformation matrix, [A ], occupies a 
central position in the theory of symmetrical 
components. It is easy to show that the matrix, 
[A ], has the following important properties 


[A]Y=[A], (167) 


where [A ] is the transposed matrix of [A ]. 
ao 
[Ay '=-[A], (168) 
n 


where [A ]~’ is the inverse of [A] and [A] is the 
conjugate of [A ]. Let us now premultiply Eqs. 
(121) and (122) by [A }-'. We thus obtain 


PIA} '(V),—pLAT MV): 
=[A}"[Z][AJAT'D. (169) 

PLA} 'U),— PLAS): 
=[A}'LYJ[LAJAT'V). (170) 
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Let us now define the matrices 


[A }(V)o=(V) 9s=(Viosr), (171) 
[AJ (Do= (Doe = (Losr), (172) 
[AT (V):=(V) e= (Vier), (173) 
[AT Di=(D ts= (eer), (174) 
[A}'(Z)](4]=([2],., (175) 
[A}°LY][A]=[Y1.. (176) 


With this notation, Eqs. (169) and (170) are 
transformed to 


P(V)os—P(V) ts =[Z](D) ts, 
PD gs— P(D ts =CY IC V) te- 


Now the whole purpose of this transformation 
is that if the matrices [Z] and [Y], have the 
ring-symmetry discussed above, then the matrices 
[Z], and [ Y], defined by (175) and (176) are 
diagonal matrices. It is easy to show that the 
matrices [Z ], and [ Y ], have the form 


(2}.=[A}(Z][4] 


(177) 
(178) 


Zi 0 00 --- O 
=| 8 8 8D car9) 
0 0 O 0 Zn 
and 
(Y}=(AS'LY ILA] 
i we ® 
ae re 
. & wa 
The elements Z, and Y, are given by 
Zr41=(ZutaZita-"Z 34+ -:- 
| +a-“—-'Z,,], (181) 
Yr41=L¥uta Vieta-*Vist+::- 
+a-@-rY,, ], (182) 


and the relations between the impedances Z;; and 
the admittances Y;; because of the ring-sym- 
metry must be substituted into (181) and (182). 
Now because of the diagonal nature of the 
matrices [Z], and [Y],, the matrix equations 
(177) and (178) separate into the m scalar 
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equations 
DP Vosr— PV ter =Zrl ter (y=1, 2, dle nN), (183) 
PI ger — Pl ter = FeV cor (184) 


Eliminating V,,, from these equations, we obtain: 


T tor = p*T gsr (p?— Y,Z,) —PY,V oer (p?— Y,2+) 


(r=1,2,--- mn). (185) 
Eliminating J;,,, we have 
V tor = P? V oor ‘(p?- Y,Z,) 
—pZ,1Iqs:/(p?-—Y,Z,-). (186) 
If we let 
Y,Z,=4a," (187) 


and make use of the transforms (135) and (136), 
we obtain 


I 4,=I sr cosh (a,2) 

— V,., sinh (a,z)(Y,/Z,)', (188) 
Ver= Vosr cosh (a,2) 

—I,., sinh (a,z)(Z,/Y,)*. (189) 


The set of Eqs. (188) and (189) may be 
written in the following matrix form 


(I),=[cosh (az) Ja(1) os 
—[sinh (az)(Y/Z)* ]a(V) os, 
(V),=[cosh (az) ]a(V) 9 
~[sinh (az)Z/V]a(Dos, (191) 


(190) 


where the matrices with the subscript, d, are 
diagonal matrices with elements given by the 
Eqs. (188) and (189) and the subscript, r, 
running from 1 to m. Now in view of Eqs. (171) 
to (176), if we premultiply the Eqs. (190) and 
(191) by [A ] and substitute from these equations 
for the matrices (J),, and (V),., we obtain 


(I) =[A ][cosh (az) JaLA FD), 


~[A [sinh (az) ¥/Z]uLA}(V),, (192) 
(V)=[A ]Lcosh (az) ][A}'(V), 
—~[A [sinh (az)Z/Y]LA}"*(),. (193) 


Equations (192) and (193) give explicitly the 
complex potential and current amplitudes in 
terms of the sending end conditions of the system 
for the case of ring-symmetry. 
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(c) Identical conductors, perfect trans position.— 
In the practical case where all the conductors are 
identical and are perfectly transposed with re- 
spect to each other and the ground plane, the 
impedance and admittance matrices of the 
system have the property that: 


Zr=Ly and Z,,=Z, 1r#s, (194) 
Yr=Y, and Y,=Y, rxs, (195) 


that is, all the self-impedance and admittance 
coefficients are identical and all the mutual 
impedance and admittance coefficients are iden- 
tical. Practically, this is a very important case 
and for this case the Eqs. (181) and (182) reduce 


to 


Z:=[(Zo+(n—1)Z], Z,=(Zo-—Z), 


r=2,3,---m, (196) 
Yi=[Yot(m—1)Y], Y-=(Yo-—Y), 
r=2,3,---m. (197) 


In this case we have only two distinct propa- 
gation constants given by the Eq. (187). This 
fact enormously simplifies the solution. 


4. Solution in terms of receiving-end conditions 


The general solutions (137) and (143) are 
formulated in terms of the sending-end currents 
and potentials. In many practical cases, the 
receiving-end potentials and currents are specified 
and the sending-end potentials and currents are 
to be determined. To treat this case, it is simpler 
to measure the distance, z, from the receiving 
end. If we follow this procedure, the funda- 
mental equations (110) and (111) remain un- 
changed in form with the exception that the left 
members are now +d(I)/dz and +d(V)/dz, 
respectively. If we make these changes in the 
equations and proceed as in Section 2 we may 
solve these equations subject to the boundary 
conditions, that (V)=(V), and (J) =(J), at z=0. 
The solutions take the form: 


(1) = [F(m,) (1), cosh (x/m,2) /D’(m,) 
r=1 


+5 [F(m,) YM), 
Xsinh (4/m,z)/D'(m,)/m, (198) 
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and 


(V) aS [G(m,) } cosh (z4/m,) /D’(m,) 
r=1 


+¥ [G(m,) ][Z]), 
r=1 
Xsinh (24/m,)/D'(m,)\/m,, (199) 


where all the quantities have the same signifi- 
cance as in Section 2, but now z is measured from 
the receiving end. By the use of these equations, 
the sending-end potentials and currents may be 
determined if the required receiving-end po- 
tentials and currents and the electrical charac- 
teristics of the line are given. The general 
solutions, (137) and (143) in terms of the sending- 
end quantities and the solutions (198) and (199) 
in terms of the receiving-end quantities may be 
put into a very convenient form in the following 
manner. Let us introduce the matrices 


[A ()]=5 [ F(m,) | cosh (z\/m,)/D’(m,), (200) 
r=1 


[B(2)]=> [F(m,) LY] 

~ Xsinh (24/m,)/D'(m,)\/m,, (201) 
[C(e)]=¥ [G(m,)] 

a Xcosh (2(m,)/D’(m,), (202) 
[D(e)]=¥ [G(m,) 2] 

~ Xsinh (24/m,)/D’(m,)\/m,. (203) 


With this notation, the general equations 
governing the distribution of potential and 
current may be written in the form 


(1) =[A (2) ](),—LB) (VY), 
(V)=(C(z) }(V),-[D(2) Do 


with z measured from the sending-end. The Eqs. 
(198) and (199) become 


(1) =[A (2) ](D), +B) (VY), 
(V)=(C(z) }(V),-+[D(2) J), 


with z measured from the receiving-end. 


(204) 


(205) 


5. Solution in terms of terminal impedances 


Equations (204) and (205) are useful in de- 
termining the current and potential distribution 
of the multiconductor line provided that the 
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sending-end or the receiving-end currents and 
potentials are specified. In some cases, these 
quantities are not specified but only the sending- 
end potentials are specified and the impedances 
of the various branches of the terminal apparatus 
at the receiving-end are given. This is shown 
schematically in Fig. 4. In this case, the receiving- 
end potentials and currents are related to each 
other by: 


(V),=(Z],-(1)-. (206) 


[Z |, is a square matrix of the mth order whose 
elements are the various self and mutual imped- 
ance operators of the terminal apparatus. To 
determine the current and potential distribution 
in this case, let us place z=s in the Eggs. 
(204), we then have 


(I),=[A(s) }(D),—-[B(s) (V),, 
(V)-=CC(s) (V)o—LD(s) Do- 


We now premultiply the first of the Eqs. (207) by 
[Z |, and equate the resulting expression to (V),. 
We then solve for (J), and obtain: 


(1),=({Z,LA (s) J4+(LD(s) ])-“(LC(s) J 
+[Z],[B(s)]})(V),. 


This determines the sending-end currents in 
terms of the given sending-end potentials and the 
known constants of the line and the terminal 
impedance matrix. To determine the potential 
and current distribution along the lines, it is only 
necessary to substitute the value of (J), given by 
(208) into the Eqs. (204). 


(208) 


VIL. CONCLUSION 


Starting with the basic electrodynamical equa- 
tions of Maxwell it has been demonstrated that 


the form of the classical differential equations of 
transmission theory are rigorously valid only in 
the case of a system consisting of perfect con- 


ductors embedded in a perfect dielectric. That is, 


the system may be specified by its self and mutual 
series impedances and shunt admittances only in 
this ideal case. In this case the propagation 


constant is given by jw/v and we have unat- 


tenuated transmission with the velocity of light 
in the dielectric medium. The wave in this case 
is a pure plane guided wave and the fields may 
be derived from one potential function. It was 
demonstrated also that for the slightly dissipative 
case of good conductors embedded in good 
dielectrics, the problem could be regarded as a 
perturbation of the ideal case. In this case, it was 
demonstrated that the system may also be 
specified by self and mutual series impedances 
and self and mutual shunt admittances and that 
the problem may be formulated in terms of these 
parameters. It was also demonstrated that these 
parameters may be calculated from the solution 
of two-dimensional potential problems. However, 
in this case the solution of the classical differ- 
ential equations leads to an equation of the mth 
order for a?, where a is the propagation constant. 
This indicates that in the dissipative case, there 
are n modes of propagation. 

Having formulated the problem in terms of 
the “‘generalized telegraphist’s’’ equation of the 
classical theory, it was demonstrated that the 
solution of these equations may be effected by the 
use of the Laplacian transformation and matrix 
The 


symmetrical cases of cables and transposed aerial 


notation in a very convenient manner. 


lines were considered and solved in general. 
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An Operational Treatment of Electromagnetic Waves Along Wires 
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Starting with the fundamental field equations of Maxwell, the equations governing the 
propagation of electromagnetic waves along two parallel wires are derived. The treatment is 
greatly facilitated by the use of the Laplacian transformation or modern rigorous operational 
calculus. First the case of waves traveling along perfectly conducting wires is considered and 
the dissipative case is treated as a perturbation of the ideal perfectly conducting one. An 
expression for the deformation of the front of an initially rectangular wave due to the “skin- 
effect”’ is derived and the manner in which the wave progresses along the wire is studied. 


INTRODUCTION 


HE literature discussing the propagation of 

electromagnetic waves along a pair of two 
parallel wires is quite extensive.'~* The type of 
analysis found in the literature may be divided 
into two classes. Some investigators attack the 
problem in a fundamental manner beginning 
with the field equations of Maxwell and concern 
themselves with studying the space propagation 
of alternating potentials and currents along the 
wires. By determining the spatial behavior of a 
simple sinusoidal wave with respect to time, it 
is then possible to study the effect of arbitrary 
functions of time and transient effects by means 
of the Fourier integral method. This method 
takes into account the variation of the param- 
eters of the equations with frequency. 

Another group of investigators base their 
analysis on the hypothesis that the two-wire 
arrangement may be specified by the generaliza- 
tion of the circuit theory concept and consider 
the system as being made up of a continuous 
chain of lumped circuits having certain values of 
capacitance, inductance, resistance, and leakage 
conductance. This type of analysis is much 
simpler than the one based on field theory and 
may be shown to be rigorously valid only for the 
case of perfect conductors embedded in perfect 
dielectrics but is a very good approximation to 
the actual case provided that the conductivity of 
the conductors is much greater than that of the 
dielectric. 

The analysis presented in this paper begins 
with Maxwell’s equations. These equations are 
transformed by the Laplacian transformation 
and the dissipationless case is solved rigorously. 


800 


The case of finite conductivity is then considered 
and solved as a perturbation of the ideal case. 
The usual ‘operational telegraphist’s equations”’ 
are then derived and the attenuation and dis- 
tortion of a wave traveling along the line is 


studied. 


I. STATEMENT OF THE PROBLEM 


Let us consider two straight cylindrical con- 
ductors parallel to one another whose lengths are 
very great compared to their distance apart. 
The two parallel conductors will be taken parallel 
to the z axis of a Cartesian reference frame as 
shown in Fig. 1, and the dimensions are indicated. 

The problem to be considered in this discussion 
is that of determining the general nature of an 
electromagnetic wave which is initiated by the 
application of a sudden impulse of potential at 
one end of this line, (s=0) at t=0. Since the 
interest is centered on determining the velocity, 
attenuation, and distortion of the wave as it 


YI 











{ 
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progresses along the parallel conductors, the 
perturbing effect of the terminals of the system 
will be neglected in the discussion. 


II. THE GENERAL FIELD EQUATIONS 
AND THE POTENTIALS 
The discussion will begin with the Maxwell 
field equations. In a rationalized m.k.s. system 
of units these equations take the form 


Vx E=-—90B/at, (1) 
VxH=1a+0D/dat, (2) 
v-B=0, (3) 
V-D=p, (4) 
B=.,H, (5) 
D=kE, (6) 
i=ck, (7) 


where, in this system of units, 

E =the electric intensity (volts/meter), 

H =the magnetic intensity (ampere turns/meter), 

B =the magnetic flux density (webers/meter?), 

D =the electric displacement (coulombs/meter?), 

o =the conductivity (1/ohm meter), 

k, =the inductive capacities of the medium, 

k,=k/ko the dielectric constant or specific inductive 
capacity of the medium, 

ur=pu/uo=the permeability of the medium, 

ko = (1/367) X 10~° farad /meter, 

uo=4X 10-7 henry /meter, 

p=the charge density (coulombs/meter’). 


The factors uo and ko satisfy the fundamental 
relation: 


C= 1(poko)' = 2.998 X 108 meters/second 
(velocity of light in vacuum). (8) 
The potentials 


The analysis of the electromagnetic field is 
facilitated by the use of auxiliary functions 
known as potentials. This may be done by 
placing 


B=VxA, (9) 
E=—-—Vo-—0A/at, (10) 
where 
@=the scalar potential (volts), 


A=the vector potential (joules/ampere). 
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In a medium of inductive capacities, k and u, 
and conductivity o the potentials satisfy the wave 
equations 


V°o — kud*/ dt? — pod /dt =0, 
V?A —kyd?A /dt?—podA /dt =0. 


(11) 
(12) 


The potentials are not independent but satisfy 
the relation: 


V-A+pkdo/dt+poo=0. (13) 
III. THe LAPLACIAN TRANSFORMATION OF THE 
WAVE EQUATIONS 


The solution of the problem under considera- 
tion will be effected provided that we find solu- 
tions of the wave equations (11) and (12) that 
satisfy the proper boundary conditions at the 
surfaces of the conductors. Since the conductors 
are taken to be in the z direction, we may assume 
that the magnetic vector potential, A, has only a 
z component, A. This assumption is rigorously 
fulfilled in the case that the medium surrounding 
the two wires is a perfect dielectric since, in this 
case, the conduction currents will have only a z 
component. In an actual case of imperfect di- 
electrics, the leakage conductance currents be- 
tween the conductors is vanishingly small com- 
pared to the current in the conductors, so that 
the above assumption is justified in all practical 
cases. We may, therefore, replace the magnetic 
vector potential, A, in the above equations by its 
z component which we write as A _ without 
the bar. 

To facilitate the analysis, we will make use of 
the methods of the operational calculus in the 
modern rigorous Laplacian transformation form 
(5). To do this, let us introduce the following 
functions: 


emotes) =Pf erated (14) 

eAlxy.p)=p J | (15) 
0 

e“Bixy.p)=p [rb ey.sidt (16) 
0 

e~* B,(x,y,p) =p fe” Boxysindt (17) 
0 

801 








Using the symbolic notation of the operational 
calculus, the above relations may be written in 


the concise form 
e-"*,(x,y,p) = o(x,y,2,t), ete. (18) 


In the above equations, a, is the “operational 


propagation constant,’’ whose value is to be 


determined. Substituting (14) and (15) into the 
wave equations (11) and (12) and making a few 
reductions, we obtain 


V*b. = (p*, ve + pou—a*)d: (19) 
anc 

V°A.= (p*, ve? + pou—a*)Az, (20) 
where vo =1/(uk)! is the velocity of light in the 


medium under consideration. The same process 
applied to the Eq. (13) gives the relation 


A,=(p/av?+ou/a) dr. (21) 


The transforms of the fields, satisfy the 
equations 


B, =i0A ,/dy—jdA,/dx, (22) 


E,.= —idd,/ dx —j0d,/dy+k(agd,— pA.), (23) 


where i, j, and k are the unit vectors along the 
x, y, and z axes, respectively. 


IV. IpeAL CASE, PERFECT CoNbDuUCTORS EM- 
BEDDED IN A PERFECT DIELECTRIC 


Before considering the general dissipative case, 
let us consider the ideal dissipationless case in 
which the two parallel conductors have infinite 
conductivity and are embedded in a homo- 
geneous dielectric medium of zero conductivity. 
In this case, since the dielectric is assumed to be 
a perfect insulator, there cannot be an x or y 
component of conducting current and it follows 
that the assumption that the magnetic vector 
potential has only a z component is rigorously 
valid. We now consider the form of the equations 
in the conducting and dielectric regions. 


1. Dielectric region 


In the perfect dielectric medium surrounding 
the two wires we have 


o=0. (24) 
Hence in the dielectric region, the functions ¢, 
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and A, satisfy the equations 


V*b.=(p* Vo — a") di, (25) 
V?A = (p?/v0? —a*) A). (26) 


Now in the dissipationless case under con- 
sideration, the energy flow must be in the ¢ 
direction since no energy is being consumed in 
the perfect conductors or in the perfect dielectric. 
It thus follows that the Poynting’s vector, 
N=(ExAj), (watts/meter?) must have only a z 
component. This requires that the transform of 
the electric intensity, E,, must not have a 2 
component. From (23) this leads to: 


ad:= pA:. (27) 


However, Eq. (21) for this case gives the re- 
lation: 


A.= pd: /av*, (28) 


where v is the velocity of light in the dielectric 
medium under consideration. Substituting this 
into (27) gives 

a= +p, Vo. (29) 


Placing this value of a in (25) and (26), we see 
that both the transforms of the potentials satisfy 
the two-dimensional Laplace’s equations, 


V°*o.=0, (30) 
V?A.=0, (31) 


in the dielectric region. In view of (29), the re- 
lation (28) takes the form 


A:=¢1, Vo- (32) 


Hence the electric and magnetic fields may be 
derived from the single potential function, ¢,, in 
the form 


— 
B.= - (idg, Oy — jog, Ox), (33) 


vo 


EB, = —id¢,/dx —jd¢,/dy. (34) 
From this we obtain the relation 
E,-B,=0. (35) 


This shows that in the case of perfect con- 
ductors embedded in perfect dielectrics, the 
electric and magnetic vectors are both in planes 
normal to the z axis and normal to each other. 
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2. The conducting region 


In the perfectly conducting region, c—* and 
the equations (19) and (20) reduce to 


V*or= poudr, 
V2A ‘= popA t- 


(36) 
(37) 


These equations are of the form of the well- 
known heat diffusion equation. Since in this 
case we are assuming infinite conductivity, the 
solution of these equations shows that there is 
no penetration of electric or magnetic field into 
the conductors. Accordingly, we may take ¢, and 
A, as constant over the conductor cross sections. 
It is thus seen that in this ideal case there are 
no fields inside the perfect conductors and that 
the fields outside the conductors are derivable 
from a single scalar potential function, ¢:, that 
satisfies Laplace’s equation (30) in the dielectric 
region and reduces to prescribed constant values 
at the surfaces of the conductors. 


3. The wave velocity 


In order to study the manner in which the 
fields and potentials are propagated along the 
ideal conducting system, it will be necessary to 
make use of some well-known theorems of the 
Laplacian transformation. In the literature of the 
operational calculus it is customary to use the 
notation 


eb 4(x,y,p) = O(X,9,540) = 


when two functions are related to one another 
by equations of the type (14). We then say that 
e-“,(x,y,p) is the direct Laplacian transform of 
$(x,y,z,t) and conversely, $(x,y,z,t) is the inverse 
Laplacian transform of e~¢,(x,y,p). Let us now 
assume that at z=0, we have the relation 


o(x,y,0,t) = do(x,y,2). (39) 


Placing z=0 in (38), we see that ¢,(x,y,p) is the 
transform of the function ¢0(x,y,t). Now from a 
basic theorem in the operational calculus, we 
have the result that if 


g(p)=h(t), (40) 
then 
( 0 t<k 
e*Pg(p)= « k, positive. (41) 
h(t—k) t>k 


Since in the ideal, dissipationless case under 
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consideration, we have a= p/v, substituting this 
relation into (38) we obtain 


. | 0 t<z, ‘Vo 
$(x,y,2,t) = 4 / (42) 
\bo(x,y,t—2/v9) t>2z/v0. 


It is thus seen that the potentials and fields are 
propagated along the conductors without attenu- 
ation or distortion with the velocity v, which is 
the velocity of light in the dielectric medium 
between the conductors. 


4. The circuit formulation of the ideal case 


(a) The coefficients of capacitance.—Having de- 
duced the general properties of the non-dissipa- 
tive case from the basic field equations, let us 
now turn to a formulation of the problem in 
terms of the circuit concept involving the coeffi- 





Fic. 2. 


cients of capacitance and inductance of the two 
conductors. In order to obtain the ordinary 
circuit formulation of the problem, it is necessary 
to write the equation of conservation of elec- 
tricity for the case under consideration. To do 
this, let us take the divergence of both sides of 
Eq. (2) and obtain 


V-(VX A) =0=V-(a+aD/a0). (43) 


Let 7;(z,t) and 72(z,t) be the currents flowing in 
conductors 1 and 2 as shown in Fig. 1 and let 
qi(z,t) and q2(z,t) be the charges per unit length 
at the surfaces of conductors 1 and 2 and let us 
introduce the transforms, 


e~*],(p) +1,(z,t), 
e~"0,(p) =4q,(z,). 


(44) 
(45) 


Then if we apply Gauss’ theorem to the 
surface bounding a length dz of the rth cylindrical 
conductor as shown in Fig. 2 of preceding paper, 
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we obtain in view of (43), 


al,= — 00 § (86x ‘On)ds, r=1,2, (46) 
where ¢;, is the value of the potential transform, 
¢o:(x,y,p) evaluated at the surface of the rth 
conductor and ga» is the operational conductivity 
of the dielectric medium defined by: 


oo=oat pha. (47) 


oa is the conductivity and k, is the inductive 
capacity of the dielectric medium surrounding 
the conductors. By an application of Gauss’ 
theorem to Eq. (4) over the same surface of 
Fig. 2 of preceding paper, we have 


Q,(p) = —ke $ (861 ‘On)ds. r=1,2 (48) 


In the ideal case under consideration, o,=0, 
and hence from (46) and (48), we obtain 


al,=pQ,. 


Now the function ¢,(x,y,p) satisfies Laplace’s 
equation (30), and reduces to the values ¢,, and 
¢di2 at the surface of the conductors No. 1, and 
No. 2, respectively, and also satisfies the relation 
(48). Hence by the ordinary methods of potential 
theory, we may express Q; and Qz in the form 


r=1,2 (49) 


O1=Cruduatcidre, 


(50) 
Q2=CoidiitCoodre, 


where the c,, quantities are the coefficients of 
capacity and induction per unit length of the 
conductors. Since conductor No. 2, serves as a 


return for conductor No. 1, we have 
I,=-I,=—-I. (51) 


Combining Eq. (49) with (50) and using the 
relation (51), we obtain: 


al = p(cidiutciedie), 


(52) 
—al = p(C2191+€22912). 
From this we have 
al =p(di— 12) (C11 — C12”) / (C1 +2e2+2¢12) (53) 
if we let 
(d1—¢2) = V(p), (54) 
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where V(p) is the transform of the potential 
difference between the conductors and related to 
the actual potential difference, v(z,t) by the 
equation 


e~* V(p)=v(z,f). 


The quantity 


(55) 


C= (€11€22 — C19") / (C1 +622 +2012) (56) 


is the coefficient of capacitance per unit length 
of the conductors. For the arrangement of two 
conductors as shown in Fig. 1, this coefficient is 
given by 


C=-rk./cosh (d/2r) (farads/meter), (57) 


where k, is the specific inductive capacity of the 
dielectric medium in which the conductors are 
embedded. If the distance separating the con- 
ductors is great in comparison to the radii of 
the conductors so that the relation d>r is 
satisfied, then (57) reduces to the approximate 
expression 


C=rk,/log, (d/r) (farads/meter). (58) 


It follows, therefore that the Eq. (53) may be 
written in the form 


al=pCV. (59) 


(b) The coefficient of inductance.—To obtain 
the relation involving the coefficient of induct- 
ance, we make use of the Eqs. (27) and (32) and 
write 


ag: = pd/Vo. (60) 


The transforms of the potential function at the 
surfaces of the conductors may be expressed as a 
linear combination of the transforms of the 
charges per unit length of the conductors. In 
view of (60), we may write 


p 
Agu =—(piQit pie), 


(61) 
p 
Adi2=—(p2101+ p22Q2), 


Vo 
where the p,, quantities are the well-known 
coefficients of potential. Now since conductor 
No. 2, is the return of conductor No. 1, it 
follows from (49), that 


Q2.= —Qi=—-Q. (62) 
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Substituting this relation into (61) and sub- 
tracting the second equation from the first, we 
obtain 
Pp 
a(o11— biz) =-Q( Pit po2— 2pi2). (63) 


Uv 


Making use of the well-known relation of poten- 
tial theory, 


(Put po2—2pi2) =1/C, (64) 
we may write (63) in the form 
aV=pQ/mC= pl/vePC. (65) 
The quantity 
Lo=1/v°?C (henry /meter) (66) 


is defined as the coefficient of external inductance 
per unit length of the system of conductors 
under consideration. In terms of this coefficient, 
we may write (65) in the form 


aV=pLol. (67) 


Equations (59) and (67) are the well-known 
equations on which the classical or circuit theory 
of wave propagation is based. It is seen that 
these equations are rigorously valid for the 
dissipationless case under consideration. 


V. Tue CASE OF FINITE CONDUCTIVITY 


Having considered the ideal, non-dissipative 
case in detail, it will now be used as a model for 
the analysis of the actual case in which the con- 
ductivity of the conductors although great is not 
infinite and that of the dielectrics although 
small is not zero. This is the condition that is 
satisfied in any efficient transmission system. 
In order to proceed with the analysis of this 
case, it is necessary to make certain approxima- 
tions which will be enumerated. 


1. General approximations 


a. The conductivity, o, of the conductors will 
be assumed to be of the same order of magnitude 
as that of copper (5X10? mhos/meter) and the 
relative dielectric constant, k,, of all the materials 
will be taken as being of the order of unity. 

b. The conductivity of the insulating medium 
will be taken to be of the order of 10-* to 10~'° 
mho/meter. 
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c. The relative permeability, y,, will be as- 
sumed to have the order of magnitude of unity 
in all dielectrics and non-magnetic materials. 

d. Because of the low conductivity of the 
dielectric medium, the longitudinal or z com- 
ponent of the magnetic field may be neglected in 
comparison with the transverse components. 
Hence the magnetic vector potential, A, will be 
assumed to have only a z component. 

e. The separation between the conductors will 
be assumed to be so great in comparison with 
the radii of the conductors that the “‘proximity 
effect”’ may be neglected. That is, the assumption 
will be made that the effect of the external field 
on the distribution of current density and the 
field within the conductors may be neglected. 

f. In view of the above approximations and 
the assumed deviation from the ideal case, it 
will be assumed that the operational propagation 
constant is of the order of p/v as in the ideal 
case. It will be the object of the analysis to obtain 
the perturbation terms of the operational propa- 
gation constant introduced by the presence of 
dissipation in the system. 


2. The general equations 


We now turn to a consideration of the Eqs. 
(19) and (20) satisfied by the transforms of the 
scalar and vector potentials and consider the 
form these equations take in the dielectric and 
conducting regions. 

(a) Dielectric region.—In the dielectric region, 
we have 


b:= (p?/v0? + poaa—a*) du, (68) 
A,= (p?/ve + poaa—a)Az, (69) 


where og and ya are the conductivity and 
permeability of the dielectric region. The assump- 
tion will now be made that o, is very small and 
that a? may be taken approximately equal to 
p*/v0?. To this degree of approximation we may 
place the right members of (68) and (69) equal 
to zero. It is easy to show from a consideration of 
the equation V-E=0 that is satisfied by the 
electric intensity in the dielectric region that the 
above assumption is consistent with the assump- 
tion that: 


(the electric transverse field in the dielectric 
region)> (longitudinal field in the same region). 
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This is also consistent with the consideration 
that the energy flow as expressed by the Poynt- 
ing’s vector, N=E XA, shall have its longi- 
tudinal component much larger than its trans- 
verse component in the dielectric region. 

(b) The conducting region.—As a consequence 
of the assumption that //,=0, it follows from the 
Maxwell equation (1) that in the conductors, the 
electric transverse field is of the order of magni- 
tude of 1/voux compared to the longitudinal 
electric field, E,. Since 


z 


this is a very small 


quantity, we are led to the condition: 


(electric transverse field in the conductors) 
«(longitudinal field in the conductors). 


From an energy consideration, this is con- 
sistent with the assumption that the dissipation 
of energy in the conductors due to longitudinal 
currents is very large compared to charging 
currents. From this it follows that 


Vo.—0 (70) 


the conductor cross sections and 
hence the function, ¢, assumes the constant 
value ¢,, throughout the cross section of the rth 
conductor. In view of the relation (23) it is 
seen that the transform of the z component of 
the electric intensity, E,, satisfies the equation 


throughout 


VE = (p?/v?+ pou—a*) En. (71) 


In view of the fact that o is a large quantity in 
the conducting region and that a? and p*/v? are 
quantities of comparable magnitude, we may 
write 


VE = popki. (72) 


This is the well-known heat diffusion or ‘‘skin- 
effect” equation. The quantities o and yu in (72) 
are the conductivity and permeability of the 
conducting medium. 


3. The circuit formulation of the dissipative case 

Having established the equations satisfied by 
the potentials and fields in the various regions, 
the circuit formulation of the problem will be 
undertaken and the operational propagation 
constant will be obtained. The procedure follows 
the same reasoning as in the ideal, non-dissipa- 
tive case. To obtain the relation analogous to 
(59) of the ideal case, we combine (46) and (48) 
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and write: 


al,=(ou4, kat p)Q-(p), r=1,2, (73) 


where o, is the conductivity and ky the specific 
inductive capacity of the dielectric. In view of 
the fact that V*¢,=0 in the dielectric and reduces 
to a constant value at the surface of the con- 
ductors as in the ideal case, we again make use 
of the relations (50) and write (73) in the form: 


al,=(¢a Rat p)(Cridiatcredre), r=1,2. (74) 
Now since J;= —J.=IT] we obtain 
al=(G+ pC) V, (75) 


where C is the capacitance per unit length of 
the conducting system as given by (57) and G is 
the leakage conductance per unit length de- 
fined by 


G=Cou Ra. (76) 


The relation (76) is valid provided the con- 
ductors are embedded in a homogeneous di- 
electric medium. For aerial lines where the chief 
contribution to the leakage coefficient comes 
from the imperfect insulation of the supports, 
the Eq. (76) is not valid and the leakage coeffi- 
cient must be determined by measurement. The 
determination of the relation that corresponds to 
(67) of the ideal case is a somewhat more difficult 
matter. To determine this relation, we make use 
of Eq. (23) expressing the z component of the 
transform of the electric intensity in terms of the 
potential transforms. This may be written in the 
form 


ago:= E,.+ pA. (77) 


Now at the surfaces of conductors No. 1 and 
No. 2, we have 
don =EiitpAn, ™ 
(78) 
Aoi. = Ep.2+ pA 2, 


where the subscripts indicate that the quantities 
are to be evaluated at the surfaces of the re- 
spective conductors. If we subtract the second 
equation from the first, we obtain 


(ou— O12) = (Etri— Erez) + p(Au—Avw). (79) 


Now by comparison with Eq. (63) of the dissipa- 
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tionless case, we may write 


a V= (Ein—Eu2)+pLol, (80) 


where Ly is the coefficient of external inductance 
defined by (66). It remains now to obtain an 
expression for the transforms of the z components 
of the electric intensities at the surfaces of the 
conductors in terms of the currents and geo- 
metrical dimensions of the system. To do this, 
Eq. (72) must be solved subject to the proper 
boundary conditions. 

The solution of the operational skin-effect equa- 
tion.—If it is assumed that the conductor separa- 
tion is so great that the effect of the external 
field on the distribution of current density and 
the fields within the conductor may be neglected, 
then the “skin-effect’’ equation (72) may be 
written in the polar form: 


1 
@E,,/dp?+-dE,./dp = popE:., (81) 
p 


where p and @ are polar coordinates taken from 
the center of the circular cross section of the 
conductor as shown in Fig. 2. The solution of 
Eq. (81) valid for a solid conductor, is 


Eu=BJ(np), (82) 


where B is an arbitrary constant, and Jo is the 
Bessel function of the first kind and zeroth 
order, and 


j=V-1. 


The transform of the total current flowing 
across the conductor cross section in terms of 
the current density, u, and the electric intensity 
transform is given by 


1= [ f udxdy= f olrpE.dp 
0 


=2roB f Ju(np)pde=2xoBrJ(nr)/n, (84) 


n= j(pou)! (83) 


where J; is the Bessel function of the first kind 
and first order. From (82), the constant, B, may 
be obtained in terms of the value of E,, at the 
surface of the conductor 1, Ey: by placing 


p=r. We thus obtain 
B = Ein/Jo(nr). (85) 
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Substituting this into (84) and solving for E;,1, 
we obtain 


Eui=nJo(nr)I/2norJ (nr). (86) 


Now if we let Ro be the resistance in ohms/meter 
of the double line we have 


Ro=2/neor’. (87) 
In terms of this quantity, we write: 
Eui=nrRoJo(nr)I/4Ji(nr)=Zi(p)I, (88) 


where Z;(p) is the internal operational impedance 
per meter of conductor No. 1. Now since we are 
assuming that conductor No. 2 is identical with 
conductor No. 1, we have 
Eu2= —Z,(p)I. (89) 
Substituting (88) and (89) into (80), we obtain 
aV=22Z,(p)I+pLol. (90) 


It is convenient to introduce the notation 


Zo(p) =2Z1(p) = (mr) Jo(nr)Ro/2Ji(mr). (91) 
We then write (90) in the form 
aV=Zo(p)I+pLol. (92) 


This is a generalization of (67) of the ideal 
dissipationless case. 
4. Attenuation and distortion of traveling waves 
Eliminating J from Eqs. (75) and (92), we 


obtain 


[a?—(Zy9+pLo)(G+pC) ]V=0. (93) 


Since for V #0, the bracket term must vanish 
and we have 
a=[((Zo+pLo)(G+pC) |'=a(p). (94) 


From (91), we obtain the relation between the 
transforms of potential and current, 


I(p) = V(G+pC)/(Zo+pLo). (95) 
Now from (55) we have the relation 
e~“ V(p)+v(z,t), (96) 


where v(z,t) is the potential difference between 
the conductors. If at z=0, a potential Eo(t) is 


807 











impressed between the conductors, and we intro- 
duce the notation 


Vol p) = Eo(t), (97) 


then, in view of (96), we have 
e~* Vo(p)=v(z,t). (98) 


If E,(t) is given, Vo(p) may be computed and 
the problem of determining the attenuation and 
distortion of the potential wave resolves itself 
into finding the inverse transform of e~**Vo(p) 
where the operational propagation constant, 
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a(p), is given by (94). The essential features of 
attenuation and distortion will be brought out 
more clearly if we assume that at ‘=0 and z=0, 
the impressed a constant electromotive force, Eo, 
is impressed at z=0. This condition may be 
realized physically by suddenly connecting a 
battery of potential Eo to the system as shown in 
Fig. 3. 

In this case, the potential variation is given 
by the relation, 


Eve~* = v(2,t). (99) 


Since Ep is a constant, we have only to compute 
the inverse transform of e~*?. 

The procedure of computing this inverse trans- 
form will be greatly facilitated by the use of two 
well-known relations called the “reciprocity rela- 
tions.’’"® These relations may be stated in the 


form: If g(p)=+h(t), then 


Limit g(p) = Limit h(t) (100) 
pon t—0 

Limit g(p) =Limit h(t), (101) 
p- toa 


provided that the various integrals and limits 
exist. As an example of the application of the 
second rule, let us determine the value that the 
internal impedance of the conducting system, 
Zo, assumes as [—>«, when a constant electro- 
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motive force has been impressed at t=0. By 
(101), this is given by 


(nr) Jo(nr)Ro 


Limit Zo(p) = Limit —— . (102) 
p60 p09 2S OJ (nr) 


Now from the properties of Bessel functions, we 
have: 
Limit Jo(x) =1, Limit Ji(x)=x/2. (103) 
Hence, since n=j(poyu)', we have from (102), 
Limit Zo(p) = Ro. (104) 


This is obvious from physical considerations 
since this is the d.c. resistance of the double line. 
To proceed with the determination of the inverse 
transform of e~**, we determine a as a function 
of p in the form 


a(p) =((Zo+pLo)(G+pC) } 
= p(LoC)'((1+G/pC)(1+Z0/pLo) ]'. (105) 
We must now obtain an explicit expression for 


Zo(p) as a function of p. To do this, let us 
write 


nr = jb (106) 
where b=r(poyu)'. We then have 
(jb) Jo(7b)Ro 


Sle mets, 
= G0) 





(107) 


We may make use of the following asymptotic 
expansions for the Bessel functions.’ 








e? 
Jo(j6) = [1+1/8b+9/12852 
(2b)! 
+75/1024b3+---+], (108) 
* rij) : [ b / 
—J (jb) = 1—3/8b—15/128b2 
j * (2nd)! | 


—105/1024b'—---]. (109) 
By division, we obtain 
Zo(p) = RoLb/2+1/44+3/16b 
+3/16b?+---+]. (110) 


Substituting this into (105), and letting (LoC)! 
=1/vo, where vo is the velocity of light in the 
dielectric medium surrounding the conductors, 
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we obtain 


p 
a(p) =—L(1+G/Cp)[1+ Ro/4Lop 
Vo 


+Rom/2Lo/p+oo(1/V/p)/pLo]]', (111) 
where 
m=r(on)! (112) 
and 
: P Ro F 
oo(1 ve=—sts 8mr/p 
+3/8m*p+---+]. (113) 
We thus obtain 
p / 
a(p) =—-[1+Rom/2Lop 
Vo 
+(Ro/4L0+G, ‘C), ‘p+GRom/2LoCp} 
+GRo/4LoCp? + $0(1// p)/pLo 
+Goo(1/V/p)/p*LoC]}*. (114) 


The radical may now be expanded in a series 
form by the use of the equation 
(1+x)§'=1+%/2—x?/8+---. (115) 


Carrying out this process, we obtain, after some 
reductions: 


a(p) =p/vot+botdo/pt+ei(1/V p), 


where 


Ros Cr} G/Lo\! 
n=—(—) (1n/2nLe)+-(—) (117) 
8 \Lo 2\C 


mfCr} 
dy=Re~(—) ’ 
4\L»o 


¢1(1/4/p) =a power series in (1/4/)). 


(116) 


(118) 


(119) 


Equation (116) gives the perturbation terms 
which are introduced in the operational propaga- 
tion constant because of dissipation. We now 
must determine the inverse transform of 


Eve“ = Eve~>*e-”? / vog—dozv Pe—2¢1_ 


(120) 


We may consider each factor of (120) in turn. 
The factor e~*°? is a damping factor with distance. 
This shows that the wave is attenuated as it 
travels along the line. The transform of the 
factor e~?7/"° is the same as that of the ideal 
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dissipationless case as given by (41). Hence the 
disturbance travels from z=0 with the velocity 
of vp as in the dissipationless case. To obtain the 
inverse transform of e~?*”?, we make use of the 
well-known relation 


e-'¥p=1—erf(k/2\/t), k>0O, (121) 


where erf (y) is the well-known “error function” 


defined by 
2 y 
erf(y)=— f exp(—u?)du, (122) 
Vm vo 


Tv 
placing k=doz in (121), we obtain 


edo P= 1 —erf(dye/24/t). (123) 


We now come to a consideration of e~**!, bearing 
in mind that in view of the fact that since the 
disturbance travels with the velocity of light in 
the dielectric medium, we are interested in the 
state of affairs for small values of t. The first 
reciprocity theorem may be used to compute the 
inverse transform of e~**! for small values of t. 
By the reciprocity theorem, the required inverse 
transform is given by: 

Limit e~*#10/Y?) = e9=1, 


p 2 


In view of the above, and of the transform (42), 
we finally obtain 


(124) 


0 t<z/v9 
v(z,t) = 


(125) 
Eve~*[1 —erf(do2z/2(t—2/v0)!) ] 


t> 2/Vo. 


The Eq. (125) shows the manner in which a 
wave of potential initiated by the application of 
a constant electromotive force Eo, is propagated 
along the conducting system. We see from this, 
that the larger bp and do, the smaller v(z,t) is at 
any given point. 


5. Variation of potential as a function of time 


The result (125) gives some interesting in- 
formation concerning the behavior of an initially 
rectangular wave as it progresses along the line. 
Let us place 


n=Vol —2, (126) 
$=dozv/00/2, (127) 
w=s//n. (128) 
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Then (125) may be written in the form: 


v(z,t) = Ege~*"[ 1 —erf(w) J, 


t>z/vo. (129) 


If we hold z constant and differentiate (129) 
partly with respect to ¢, we obtain: 


Ov 
—(z,t) = Ey exp( — doz) suo exp(—w*)/1/an!. (130) 
at 


This gives the slope of the potential function at 
a fixed point z, as a function of time. The slope 
of the potential function is a maximum for 

dv 

—(z,t) =0. 

i 


Carrying out this differentiation, we obtain 


dv Ey exp(—doz)svo? exp(— w*) 





ot? V4 


X [s?/n7/2?-—3/2n52]=0. (131) 


This vanishes for 


n = 2s? /3. (132) 


Substituting this value of m and s given by (132) 
and (127) into (126) gives the relation 


t= 2) ‘votd,?2? 6 =2 ‘Votto. ( 133) 


This shows that the maximum potential slope 
with respect to time occurs at any given point 
at a time (d,’2*/6) seconds after the wave has 
first arrived at the point. This time may be 
written in the form 


to =pC2?/24n*r’oLo. (134) 


This shows that the time it takes the potential 
slope to reach its maximum is greater for large 
values of z than for small ones and varies directly 
as the square of the surge admittance C/L» of 
the line. Substituting (132) into (128), we obtain: 


dv 
“) = Eye~*e~hyg/4/m(2/3)'s?. (135) 
ot max 


This may be written in the form 
ov 


~) = Eye-e-14/4/m(2/3)ide2s?. (136) 
ot max 
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It is thus apparent that at a given point z, the 
steepest slope of the wave front is inversely 
proportional to the square of the distance from 
the origin and that the steepness is proportional 
to Lo/C. Hence it is much smaller for a cable 
than for an aerial line. It is of interest to deter- 
mine the value of the slope of the wave at its 
instant of arrival. This may be determined by 
placing n—0 in (130). Carrying out the limiting 
process, we obtain 


ov 
=) =0 
ot t=z/v9 


That is, the initial slope of the wave is zero. 
Placing w=s/4/n=4/3/2 in (125), we obtain the 
value of the potential at the point when 6v/d¢t is 
a maximum to be 


(137) 


v= Ege~*(1 —erf(4/3/2) ] 
v= 0.083 Eye~>™. 


(138) 
(139) 


or 


In the same manner, the variation of potential 
along the line at a given time ¢, may be studied. 


VI. CONCLUSION 


The natural manner by which interesting and 
useful results in the rigorous analysis of traveling 
waves on wires may be obtained has been demon- 
strated. The principal result derived is the 
equation (125) giving the velocity of propaga- 
tion, attenuation and distortion of an initially 
rectangular wave traveling along the line. By the 
use of this result and the superposition theorem, 
the behavior of waves of a general form may be 
easily studied. 
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Temperature Dependence of Young’s Modulus and Internal Friction 
of Lucite and Karolith 


Joun S. RINEHART* 
Department of Physics, State University of Iowa, Iowa City, Iowa 


(Received August 5, 1941) 


The Young’s modulus and specific energy loss of Lucite and Karolith were measured at 


different frequencies in the neighborhood of 50 kc/sec. as a function of temperature. The 
temperature range for Lucite was from —55°C to 65°C, and for Karolith, 25°C to 110°C. It was 
found that the reciprocal of Young’s modulus versus temperature curves showed positive 
curvatures at higher temperatures for both materials. At room temperature the Young’s 
moduli of Lucite and Karolith were 4.72 and 6.1410" dynes/cm?, respectively. The specific 


loss was found to vary with frequency for both materials and varied in a linear fashion for 


Lucite at room temperature. Absorbed water seemed to have a pronounced effect on the 


Young’s modulus of Karolith. 


INTRODUCTION 


ECENTLY a considerable amount of inter- 

est has been shown in the processes which 
account for the internal frictional losses occurring 
in plastics materials when subjected to vibra- 
tional stresses.'* These processes are apparently 
closely connected with the plasticity and dielec- 
tric properties* of the material and more experi- 
mental data would be helpful in arriving at a 
satisfactory theory. The increased commercial 
use of plastics makes it of technical significance 
to know how the elastic properties of these 
materials depend upon temperature. In the 
present study, a dynamic method was used to 
study the variation of internal friction and 
Young’s modulus as a function of both frequency 
and temperature of two thermoplastic‘ plastics 
having a rather wide softening range. It is 
believed that the results obtained here will be of 
both technical and theoretical interest. 


* This work was done while the writer was a member of 
the 1941 summer staff of this institution. He is regularly 
a member of the staff of the Department of Physics, 
Wayne University, Detroit, Michigan. 

1 A. Gemant and W. Jackson, Phil. Mag. 23, 960 (1937); 
A. Kimball and D. Lovell, Phys. Rev. 30, 948 (1927); 
A. Gemant, J. App. Phys. 11, 647 (1940). 

2 A. Gemant, J. App. Phys. 10, 508 (1939). 

3A. Gemant, Trans. Faraday Soc. 31, 1582 (1935); 34 
269 (1938); J. App. Phys. 9, 730 (1938). 

*In general plastics fall into the two rather distinct 
classes: thermoplastic and thermosetting. The former 
have an extended softening range and behave similar to 
glass. The latter harden because of a chemical reaction. 
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MATERIALS INVESTIGATED 


The two plastics materials studied were Lucite, 
an acrylate plastic manufactured by E. I. du 
Pont de Nemours and Company, and Karolith, a 
white casein plastic manufactured by the Ameri- 
can Plastics Corporation.® These two were chosen 
because of their extensive use in industry and, 
since both are thermoplastic, they have a wide 
softening range, and hence are very suitable for 
studying internal losses during the softening of 
the material. 

The method used required specimens having a 
square cross section about 3 mm on a side and 
from 2 to 7 centimeters long. These specimens 
were machined from large commercial sheets of 
the material. 


MEASUREMENT OF YOUNG'S MODULUS AND 
SpeciFic Loss 


Young’s modulus of the specimens was meas- 
ured by a dynamic method which has been 
adequately described* and was previously used 
by the writer.’ Longitudinal vibrations wére set 
up in a composite cylinder made by cementing* 





5 The writer wishes to thank this company for Supplying 
him with a wide variety of samples. 

» S. L. Quimby, Phys. Rev. 25, 558 (1925); J. Zacharias, 
Phys. Rev. 44, 116 (1933); L. Balamuth, Phys. Rev. 45, 
715 (1934); F. Rose, Phys. Rev. 49, 50 (1936). 

7J.S. Rinehart, Phys. Rev. 58, 365 (1940); Phys. Rev. 
59, 308 (1941). 

& Several kinds of cement were tried and the results 

seemed to be independent of the cement used, as is to be 


811 








the specimen to a suitably cut square cross- 
sectioned quartz rod. The basis of the method is 
that, as the frequency of a sinusoidally varying 
potential difference, applied to electrodes at- 
tached to the quartz crystal, is varied in the 
neighborhood of the mechanical resonance fre- 
quency of the composite bar, the current in the 
circuit containing the quartz crystal passes 
through first a maximum and then a minimum as 
the frequency is increased. This mechanical 
resonance frequency, fo, lies between the fre- 
quency, fmax, Corresponding to a current maxi- 
mum, and the frequency, fmin, corresponding to a 
current minimum and is given by® 


Tate 
= max t ( min~ max) eo 
Se S J J Imex t+ I mia 


where Imax and Jmin are the maximum and 
minimum values of the current, respectively. It is 
customary to match the individual resonance 
frequencies of both the quartz and the specimen 
to within a few percent of fy. This was done in the 
present case. The difference between fmax and 
fmin, for the materials studied here ranged from 
about 500 to 2000 cycles per second, the current 
through the crystal circuit between 50 and 100 
microamperes, and the difference between the 
maximum and minimum values of the current 
from 20 to 2 microamperes, depending upon the 
amount of damping in the specimen. Three 
different quartz crystals were used and these had 
natural frequencies of 44,850, 60,110, and 68,380 
cycles per second, respectively. The frequency 
was measured with a piezoelectric clock'® by 
means of which absolute frequencies can be 
measured to one part in 100,000, and relative 
frequencies to within about one cycle per second. 
In general, however, plastics have a very broad 
resonance peak and the limit of accuracy was 
determined primarily by an inability accurately 
to locate this peak. The breadth of the 
_ peak caused an uncertainity in the quantity, 
(fmin —fmax), Of from 30 to 60 cycles per second. 


expected [see M. A. Durand, Phys. Rev. 50, 449 (1936)]. 
The most suitable for work above room temperature was 
found to be a cellulose acetate cement used in model air- 
plane construction. At low temperatures para-rubber dis- 
solved in Vaseline was used. 
*S. Siegel and S. L. Quimby, Phys. Rev. 49, 663 (1936). 
”S. L. Quimby, Phys. Rev. 39, 345 (1932). 
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If the resonance frequencies of the quartz 
specimen are fairly closely matched, then the 
resonance frequency of the specimen alone, f,, is 
given by" 


fe=fo—(fa—fo)(me/ms), 


where f, is the resonance frequency of the quartz 
alone; m,, the mass of the quartz; and m,, the 
mass of the specimen. If the specimen is long 
compared to its cross section—a condition well 
fulfilled in the present work—then Young’s 
modulus is related to the resonance frequency of 
the specimen alone by the formula 


E=(4/n*)pL*f,? 


where p is the density of the specimen, L its 
length, and m the number of half waves in the 
specimen. In most of the present measurements 1 
was equal to three. This number of half waves 
was used so that the mass of the specimen would 
more nearly equal that of the quartz ; however, to 
serve as sort of a check, one to four half waves 
were used. The results were in close agreement 
with one another. 

If the internal losses of the specimen material 
are large compared to those of the quartz™ the 
present method may be used to calculate the 
specific loss, 0, of the specimen.” This is given by 
the formula® 





Plan dotin { (fo—fmax)(fmin — fo) }* 
fa | 


where the symbols are the same as those used 
previously. 

To obtain Young’s modulus and the specific 
loss, then, it is only necessary to run a current- 
response against frequency curve and apply the 
preceding formulas. 


" This is an approximate solucion of the equation 


m.f,tan (=) +m,f_ tan (=) =0. 


fa fa 
(See reference 6.) 

12 See reference 9. A condition which holds here since 
the losses for plastics are many times those for quartz. 

18 Several different constants may be used to express the 
damping produced by the specimen; however, @ seems to 
be of most physical significance (see reference 2) and is 
defined as the fractional part of the total energy lost 


per cycle. 
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The runs above room temperature were made 
by placing the composite bar and its holder at the 
center of a long Alundum cored furnace. The 
temperature was measured by means of a 
carefully calibrated Chromel-Alumel thermo- 
couple placed close to the specimen material. To 
eliminate any lag between the temperature as 
read by the thermocouple and that of the 
composite bar, the furnace current was changed 
by steps and sufficient time allowed to elapse 
until a steady state was reached, after which the 
resonance curve and temperature were measured. 
Measurements made with increasing and de- 
creasing temperatures indicated no detectable 
lag. 

The measurements at low temperatures were 
made by using a cryostat consisting of a suitably 
insulated test tube into which the composite bar 
and holder could be placed and the whole 
immersed in a bath of the low temperature 





27 








Reciprocal Young’s modulus 








15°50 25 O25 50 75 
Temperature °C 


Fic. 1. Reciprocal of Young’s modulus against 
temperature for Lucite. 


mixture. The bar was hung vertically and sup- 
ported by fine threads. It was found, however, 
that the internal damping of the specimen was so 
great that no particular caution was necessary in 
the form of support used either at low or high 
temperatures. Two cooling mixtures were used: 
ice and water, and dry ice and acetone. In each 
case the specimen was allowed to remain in the 
cryostat until it seemed to have reached tempera- 
ture equilibrium. The temperature close to the 
bar was measured by means of a copper-con- 
stantan thermocouple. 

The densities of Lucite and Karolith were 
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Fic. 2. Specific loss in Lucite against temperature. 
Upper curve is for a specimen of mass, 0.459 g; lerfgth, 


4.47 cm; fo~68 kc/sec. Lower curve is for a specimen of 
mass, 0.781 g; length, 6.26 cm; fo~45 kc/sec. 


measured and found to be 1.18 and 1.33 g/cc, 
respectively. The lengths were determined using 
a comparator. 


RESULTS 
Lucite 


The results of measurements made on Lucite 
are shown in Figs, 1, 2, and 3. The reciprocal of 
Young’s modulus against temperature from 
about —50°C to about 60°C is plotted in Fig. 1. 
Above 60 degrees the damping was so great that 
measurements could no longer be carried out. 
This curve is approximately straight up to 0°C. 
Above zero degrees it shows a positive curvature 
which may indicate that even at fairly low 
temperatures, Lucite is already in a softened 
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Fic. 3. Specific loss in Lucite at different temperatures 

as a function of frequency. At the three temperatures, 


—55, 0, and 61°C, the assumption was made that the 
specific loss varied linearly with frequency. 
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TABLE I.- 


TEMPERATURE FREQUENCY SPECIFIC 
(°C) (CYCLES/SEC.) Loss 
—55 48,960 0.169 
—55 63,380 0.232t 
—55 72,980 0.202 

0 46,850 0.169 

0 61,000 0.1857 

0 69,130 0.248 
26 45,940 0.191 
26 59,830 0.255 
26 68,030 0.296 
61 44,320 0.375 
61 58,180 0.394f 
61 65,910 0.501 


+ The accuracy of these values is doubtful. See discussion in paper. 


Lucite. 


EQUIVALENT VIs- 
COSITY (POISES) 


Younc's Moputus 


% CHANGE IN np* 
(DyYNES/CM?) 


PER KILOCYCLE 


6.55 x 10" 


5.7 x 108 
—0.9TT 

4.6 
4.8 5.25 

0.0 
4.8 
5.0 4.72 
5.1 0.2 
+ 
8.3 3.88 

0.5 
7.5 


tt The minus sign indicates that the equivalent viscosity decreases with increasing frequency. 


condition. Young’s modulus was measured at 
several different frequencies and was found to be 
independent of the frequency. Points on the 
curve in Fig. 1 were obtained from measurements 
made at approximately 45,000 and 70,000 cycles 
per second. No detectable change in length was 
noted in the material after several temperature 
runs had been made. 

Figure 2 shows the dependence of the specific 
loss, at two different frequencies, upon tempera- 
ture. The lower curve is for Lucite run at about 
45,000 cycles per second, the upper for about 
70,000 cycles per second. Both curves are drawn 
from data taken on several runs. Different points 
on the same curve do not necessarily correspond 
to exactly the same frequency since the resonance 
frequency of the composite bar, fo, changes with 
temperature; however, for a given run the 
deviation of fo from its average was never greater 
than about 3 percent. This deviation is probably 
within the uncertainity of the measurement of 
the specific loss. It is evident from Fig. 2 that the 
specific loss increases rapidly with temperature 
from about 0°C on up. Below this temperature it 
remains approximately constant. The reason for 
this increase is undoubtedly due to the increasing 
plasticity of the material; but, for any exact 
understanding of the true nature of the processes 
of internal friction involved here, these data will 
probably have to be correlated with plastic flow 
data. It is seen that at the higher frequency, @ 
begins its rapid increase at a somewhat lower 
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temperature than at the lower frequency. Another 
run was made on the material at about 60,000 
cycles per second. The curve for it fell in be- 
tween the two shown in Fig. 2. 

In Fig. 3, the specific loss at four different 
temperatures has been plotted against frequency. 
At —55°C it is almost independent of frequency. 
It is seen to become more frequency dependent 
with increasing temperature. At room tempera- 
ture the loss varies almost linearly with fre- 
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Fic. 4. Reciprocal of Young’s modulus against 
temperature for Karolith. 


quency. If the curves for both room temperature 
and 61°C are extrapolated back, they go almost 
through the origin, thus indicating that @ is 
nearly proportional to the frequency in this 
range. This is in qualitative agreement with work 
at much lower frequency done by Gemant? on 
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paraffin wax but is very different from the 
behavior of metals for which the specific loss has 
been found to be independent of frequency.!* '5 

In discussing internal friction of plastics, 
Gemant? uses the symbol 7,,* to denote an equiva- 
lent viscosity which he considers to be de- 
termined by the solid friction, elasticity, and 
plasticity, together, of the material. The equiva- 
lent viscosity is given by 


np* =E0/4x*fo. 


‘With increasing temperature 7,* depends more 
and more on the plasticity and should become 
independent of the frequency. Table I gives the 
values of Young’s modulus, specific loss, and 
equivalent viscosity for Lucite for three different 
frequencies at 0°C, 25°C, and 61°C. The depend- 
ence of the equivalent viscosity on both tempera- 
ture and frequency is shown in the last column 
which shows the percent change in 7,* per 
kilocycle. These values have been calculated 
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Youngs modulus-cm*/dynexlO™' 
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Fic. 5. Effect of absorbed water in Karolith on Young’s 
modulus. The positive direction of the ordinate indicates 
decreasing mass. 


using the two extreme frequencies. Only these 
two frequencies were used since several runs 
were made at these frequencies while the material 
was run only once at 60,000 cycles per second 
both above and below room temperature, and the 
data of the two extreme frequencies are probably 
more reliable. The point corresponding to ap- 
proximately 60,000 cycles per second in Fig. 3 
was checked several times and is probably 
equally as good as the other two. The increase of 
the equivalent viscosity, noted here, will be 
discussed later in connection with similar data on 
Karolith. 


4 T. Read, Phys. Rev. 58, 371 (1940). 
1 R. L. Wegel and H. Walther, Physics 6, 141 (1935). 
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Karolith 


The reciprocal of Young’s modulus of Karolith 
is plotted against temperature in Fig. 4. Its 
behavior from room temperature up to about 
75°C is similar to that of many pure metals. 
Above this temperature the curve begins to 
curve rapidly upward. This upward trend is 


0.4 





us 


Specific loss 
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20 40 60 80 100 
Temperature °C 


120 


Fic. 6. Specific loss in Karolith against temperature. 
Upper curve is for a specimen of mass, 0.537 g; length, 
4.92 cm; fo~59 kc/sec. Lower curve is for a specimen of 
mass, 0.651 g; length, 6.66 cm; fo~45 kc/sec. 


much more pronounced than in the case of 
Lucite and is very likely due to the softening of 
the material since the published"® softening range 
of casein plastics is 66°C-79°C. The change in 
Young’s modulus with temperature for both 
Lucite and Karolith is much more pronounced 
than it is for most metals. The change for 
Karolith in a 100 degree temperature range is 
almost 50 percent. 

One difficulty encountered during these runs 
arose from the fact that casein plastics will 
readily absorb water. It was found that during an 
individual temperature run the specimen would 
lose considerable weight and at the same time a 
marked increase in Young’s modulus would 
occur. The density of the specimen was measured 
at intervals and remained appreciably constant.'” 
This observed change in Young’s modulus with 
water content is shown in Fig. 5. It was found, 
however, that the slope of the reciprocal of 
Young’s modulus against temperature curve was 


16H. Simonds, Industrial Plastics (Pitman, 1939), p. 161. 
177t was assumed that this loss of weight was due to 
the loss of absorbed water since it is known that the in- 
crease in volume of casein plastics is roughly the same as 
the volume of water absorbed. (See reference 16, page 182.) 
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TEMPERATURE FREQUENCY SPECIFIC 
(°C) (CYCLES/SEC.) Loss 
25 45,960 0.334 
25 61,510 0.384 
75 44,610 0.222 
75 59,450 0.265 
110 43,340 0.305 
110 58,240 0.384 


independent of the amount of absorbed water, 
provided that the runs were made in a sufficiently 
short period of time so that very little water was 
lost. By baking the material for several days at a 
temperature slightly higher than any attained in 
the temperature run, reproducible results could 
be obtained. The curves in Figs. 4 and 6 were 
obtained from material so treated. 

In Fig. 6, the specific loss of Karolith is plotted 
against temperature for runs made at approxi- 
mately 60,000 and 45,000 cycles per second. The 
upper curve is for the high frequency and the 
lower curve for the low frequency run. It is 
evident that the specific loss, for both frequencies, 
decreases, reaches a minimum, and then begins to 
increase rapidly with increasing temperature. It 
is interesting to note that this rapid increase 
starts at about 75°C, approximately the same 
temperature at which the reciprocal of Young’s 
modulus begins to show a rapid increase (see 
Fig. 4). It seems likely then that the sudden 
increase of the specific loss is due primarily to the 
increased plasticity of the material. In the middle 
portion of the curves, the specific loss shows the 
same general increase with frequency for Karolith 
as it does for Lucite. At the extreme left side of 
Fig. 6, the two curves seem to be nearing each 
other and might very well merge into one if 
carried to sufficiently low temperatures. At higher 
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TABLE II.—Karolith. 





EQUIVALENT VIs- 
CosITY (PoISEsS) 


Younc’s Mopu.us 


YJ CHANGE IN np* 
(DyYNEs/CM?) 


PER KILOCYCLE 


11.3108 6.14 10” 1.0 
9.7 
6.4 5.12 0.7 
5.8 
7.6 4.25 0.5 
71 
temperatures the two curves appear to be 


approximately parallel to each other. 

The effect of frequency upon the equivalent 
viscosity is shown for three different tempera- 
tures in Table II. The percent change of the 
equivalent viscosity is listed in the last column of 
the table. It is clear that the dependence of 
np* on frequency decreases with increasing 
temperature. This is just opposite to the behavior 
of Lucite in this temperature region. 

From the general behavior of the curves in 
Figs. 2 and 6 and the data in the last column of 
Tables I and II, the writer is inclined to feel that 
the behavior of Lucite in the temperature range, 
— 50°C to 0°C, may correspond somewhat to the 
behavior of Karolith in the approximate temper- 
ature range, 70°C to 80°C. It may be that if 
Karolith were run at a much higher temperature 
the equivalent viscosity would begin to increase 
in a fashion similar to that for Lucite; for, the 
measurements made on Lucite at low tempera- 
tures show its equivalent viscosity to be de- 
creasing with increasing temperatures as is the 
case for Karolith. 

In conclusion, the writer wishes to thank the 
Physics Department of the State University of 
Iowa for the facilities placed at his disposal, and 
Professor E. P. T. Tyndall for much help in the 
interpretation of the data. 
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Precipitation Hardening in the Fe-W System* 


CHARLEs S. SMITHt 
Massachusetts Institute of Technology, Cambridge, Massachusetts 


(Received June 7, 1941) 


An x-ray study has been made of precipitation-hardened, 7.16 percent tungsten, iron-tungsten 
alloys. It is concluded that the precipitation reaction is of the discontinuous type. On the basis 
of diffraction evidence a mechanism for the formation and growth of the precipitate is proposed 
and a picture of the crystallographic form of the precipitate deduced. 


I. INTRODUCTION 


ECENTLY x-ray studies have been made of 

precipitation-hardened alloys of Al-Cu! and 
Al-Ag? in which anomalous diffraction effects 
have been observed, caused by relaxation of the 
third Laue condition. The two-dimensional 
segregation of heavy solute atoms, Cu or Ag, has 
been advanced as a possible physical situation 
which would give rise to such scattering. An 
X-ray investigation was undertaken in the Fe-W 
system in order to study the possibility of similar 
effects in a system with a different solvent 
constituent. The search for such two-dimensional 
diffraction effects was completely negative but 
the x-ray evidence obtained was sufficient to 
allow a picture to be formed of the process of 
precipitation from a solid solution of W in Fe. 


II. MATERIALS AND PREPARATION 


The Fe-W constitutional diagram on the iron 
rich side is shown in Fig. 1. The \(Fe2W) com- 
pound has been confirmed by several workers. 
There is still some doubt about the ¢ compound 
given here as Fe;W,.. The diffraction pattern of 
the « phase was studied by Arnfelt* who from 
chemical analysis decided that Fe;We was the 
best composition and that the unit cell was 
hexagonal with a=4.73A, c=25.8A and 40 atoms 
in the cell. Westgren and Arnfelt* later using the 


* Part of a thesis presented to the Massachusetts Insti- 
tute of Technology in partial fulfillment of the require- 
ments for the degree of Doctor of Science, October, 1940. 

+t Now at the Physics Department, University of Pitts- 
burgh, Pittsburgh, Pennsylvania. 

1 A. Guinier, Ann. de physique 12, 161 (1939). 

2C. S. Barrett and A. H. Geisler, J. App. Phys. 11, 733 
(1940). 

3H. Arnfelt, Carnegie Scholarship Memoirs (Iron and 
Steel Institute 17, 1 (1928). 

4H. Arnfelt and A. Westgren, Jerkontorets Annaler 
(1935), p. 185. 
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same data, decided that a better unit cell was the 
equivalent rhombohedral one with a=9.02A and 
a=30.5° having a volume 4 as great as the 
hexagonal one, and now with 13 atoms per unit 
cell. The only formula consistent with this and 
with the density is Fe;W., leading to unique 
atomic positions which were verified with inten- 
sity measurements. The metallographic work® ® 
finds no single phase at the composition Fe;W.g, 
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Fic. 1. The Fe-W constitutional diagram 
on the iron rich side. 


5 W. P. Sykes and K. R. Van Horn, Trans. A.I.M.E. 105, 
198 (1933). 

6 E. P. Chartkoff and W. P. Sykes, Trans. A.I.M.E. 89, 
566 (1930). 


817 








-— 


but does find one at Fe;W:s with some homo- 
geneity range. Apparently the situation is similar 
to that of CuAl, where the compound exists only 
in a composition range outside the stoichiometric 
formula. 

The alloy{ studied was a 7.16 atomic percent 
tungsten material prepared by melting the pure 
iron and tungsten powders together. The material 
was rolled down to a ribbon }”’ wide and 0.001” 
thick and heat treated in this form. 

Reference to the constitutional! diagram, indi- 
cates that solution heat treatment temperatures 
greater than 1380°C are necessary. Solution heat 
treatment was carried out at 1400°C in a glowbar 
furnace and hydrogen atmosphere for a minimum 
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Fic. 2. The lattice parameter of the body-centered cubic 
solid solution of Fe and W as a function of W content. The 
straight line is drawn to connect the parameter of pure Fe 
and pure W. The circles are experimental points. 


of 30 minutes. Quenching was accomplished by 
simply removing the samples to the air, the thin 
ribbons cooling to blackness within 5 seconds. 
Samples were aged for }, 3, 1, 2, 5, 10, 25, 50, 
and 100 hours at 600, 650, 700 and 725°C. The 
furnace temperature was continuously controlled 
to 2°C. Samples were placed in special high 
temperature glass tubes, which were evacuated 


with a fore pump, and sealed off in order to avoid 


oxidation during aging. 


t All the alloys used in this investigation were made and 
supplied by Mr. W. P. Sykes, Cleveland Wire Works, 
General Electric Company. The author wishes to acknowl- 
edge the assistance of Mr. Sykes who has done a large 
part of the previous work which will be referred to, on the 
constitutional diagram of iron and tungsten and on pre- 
cipitation hardening. 


818 





III. Latrice PARAMETER CHANGES 


The composition of the solid solution was 
followed throughout the aging by means of lattice 
parameter measurements made with a_ back 
reflection focusing camera and FeKa radiation. 
C.P. annealed copper powder was used to cali- 
brate each film for shrinkage. Because of the large 
grain size created at the high solution heat 
treatment temperature, it was necessary to 
abrade the sample with emery in order to obtaina 
smooth diffraction line, thus reducing the accu- 
racy of the measurements to 0.0005A. A separate 
portion of each sample was used for these 
measurements. 

It was found desirable to make a redetermi- 
nation of the lattice parameter—composition 
curve. For this purpose alloys of 0.00, 1.63, 4.90 
and 7.16 atomic percent tungsten were solution 
heat treated, quenched and measured. The 
results are in Fig. 2. 

A lattice parameter versus time curve for a 
series of aged alloys is shown in Fig. 3. The 
double ordinates for certain time values represent 
a definitely split 220 diffraction line. Other lines 
observed in a side reflection focusing camera 
showed this effect also; consequently the cause of 
the effect is interpreted as the coexistence of two 
solid solutions of different composition. This 
effect, which from the point of view of a study of 
the matrix is evidence of precipitation, starts 
after hardening has commenced and ends just as 
hardening is completed?’ for each aging tempera- 
ture studied. 

The final lattice parameter attained in each of 
the four series of specimens was essentially the 
same, independent of aging temperature, and 
corresponds to a composition of W=1.4 atomic 
percent. This is less than the published equi- 
librium value of 3.0 atomic percent, determined 
by microscopic methods. Since the present lattice 
parameter method is in fact a direct determi- 
nation of the equilibrium composition at these 
temperatures, the value of 1.4 percent is probably 
more nearly correct. 

This type of parameter-time curve is charac- 
teristic of discontinuous precipitation® in which 


7W. P. Sykes in A.S.M. Age Hardening of Metals, Sym- 
posium, 1939. 

8 R. F. Mehl and L. K. Jetter in A.S.M. Age Hardening 
of Metals, Symposium, 1939. 


JOURNAL OF APPLIED PHYSICS 











“n” 


4 


Leoftrce Parameter 





Time 


oY Aging - Hours 


Fic. 3. The lattice parameter-time curve for the material 
aged at 700°C. Other aging temperatures result in a similar 
type of curve with the same initial and final values of lattice 
parameter, but the change-over time is later the lower the 
temperature. 


the reaction 
a;—a;+ precipitate 


is supposed to go to completion in restricted 
regions and to remain at the original value in 
other regions, the reaction proceeding by the 
growth of the completely precipitated regions. 
Such restricted regions have been found in other 
cases to be in the neighborhood of either slip 
planes or grain boundaries. 

Since Sykes’ photomicrographs’ of aged ma- 
terial show clearly that in the present case these 
regions are not grain boundaries, they are there- 
fore presumed to be slip planes. In view of the 
fact that no well-defined slip planes exist in body- 
centered cubic metals and that the material 
received no intentional working previous to aging, 
the situation is not clear. Work which would 
resolve this point would be the determination of 
the crystallographic orientation relation between 
precipitate and matrix. 

The parameter-time curves show that the 
reaction goes to completion. Comparison with 
Sykes’ hardness-time curves’ shows that maxi- 
mum hardness occurs substantially at the same 
time as the completed reaction. Further, Sykes 
has observed that the greater the tungsten 
content of the original solid solution, the greater 
the maximum hardness obtained. Therefore the 
conclusion that the hardness of the aged alloy is 
due merely to the presence of a definitely 
precipitated phase seems warranted. 
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IV. PRECIPITATED PHASE 


In most precipitation hardenable alloys, x-ray 
methods have failed to detect the precipitating 
phase until over-aging has occurred. In the 
present system Sykes has detected diffraction 
lines due to a second phase at the time corre- 
sponding to maximum hardness.* Because tung- 
sten has a high atomic number relative to that of 
iron, the diffraction pattern from a tungsten rich 
compound is strong and should be readily de- 
tectable. With a refinement of technique, even a 
study of the precipitating phase should be 
possible. 

The arrangement employed is indicated sche- 
matically in Fig. 4.9 FeKa radiation is mono- 
chromated by a bent rocksalt crystal, which 
focuses the radiation to a fine line focus. If a 
transmission sample, the film, and the focus of 
the monochromator are placed on the circumfer- 
ence of a circle, focusing conditions obtain and 
sharp diffraction effects may be observed. 
Monochromation is essential for elimination of 
background due to the excitation of fluorescence 
FeKa radiation by the continuous radiation 
present in the x-ray beam. The camera was also 
evacuated to eliminate air scattering. Further- 
more, the arrangement is suitable for detecting 
diffraction at small sin @/A, which is important 
in the present case because the intermetallic 





Fic. 4. A schematic drawing of the monochromator 
camera arrangement used to study the diffraction pattern 
of the precipitating material. 


compounds concerned can best be differentiated 
at small angle where the low order diffraction 
lines from the large unit cells occur. 

The optimum thickness 1/y of an iron sample 
for FeKa radiation is 0.0004”. To attain this 
thickness the 0.001” ribbon was electrolytically 
polished with Jacquet’s solution” in order to 
avoid cold working the samples. 

® Described fully in C. S. Smith, Rev. Sci. Inst. 12, 312 


(1941). 
19 P, Jacquet, Comptes rendus 208, 1012 (1939). 
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Typical diffraction patterns resulting from this 
work are reproduced in Figs. 5 and 6, and are 
graphically indicated in Fig. 7 together with the 
diffraction patterns of the intermetallic com- 
pounds Fe,.W and Fe;W<g. In these patterns the 
ring of dense spots near the outside is the solid 
solution 110 line; all other maxima are caused 
by the precipitated phase. These other effects are 
as follows: in Fig. 5, a ring of sharp faint spots at 
d=2.37A, a very faint broad band at d= 3.94A to 
d=4.40A and a broad band at d=2.00A to 
d=2.18A;in Fig. 6 the same sharp ring at d= 2.37 
now more intense, the first-mentioned band 
sharpened to a broad line at d=4.20A and the 
second band sharpened to three broad lines at 
d=2.17, 2.10 and 2.01A. 

The precipitate diffraction pattern was first 
detected after aging for 20 hours at 600°C, 1 hour 
at 650°C, and 3 hour at 700°C, thus confirming 
the above conclusion that precipitation is con- 
comitant with hardening. Comparison of Fig. 6 
with the upper part of Fig. 7, shows that there 





Fic. 5. The diffraction pattern of the alloy aged for 5 hours 
at 725°C. 





Fic. 6. The diffraction pattern of the alloy aged for 50 
hours at 725°C. These patterns cover the range of d=50A 
to d=1.90A (@=1.5° to @=30°) with FeKa radiation. The 
ring of over-exposed spots at the edges of the pattern is 
the solid solution 110 line. 
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are definite similarities between the pattern of 
the precipitate and those of Fe;W, and Fe.W but 
that there is no basis for deciding between Fe;W, 
and FesW nor indeed for proposing any sort of 
strictly crystalline structure from the pattern. 
Since, however, there is a definite similarity 
between the precipitate pattern and that of 
Fe;W, and Fe.W, a study of the latter two 
structures has made it possible to form an idea of 
the structure of the precipitate and the mecha- 
nism of its formation. 

Figure 8 shows the crystal structures of FesW 
and Fe;W.. FeeW is hexagonal with a=4.73A, 
c=7.70A, and Fe;W, may be referred to hexago- 
nal axes with a=4.73A and c=25.8A. In these 
structures the common feature, giving rise to 
identical a parameters and therefore to common 
hko d distances, is a plane network of iron atoms 
parallel to the basal plane. These networks, 
shown in Fig. 9, occur at the heights z=1/4, 3/4 
in FeoW and at z=1/12, 3/12, 5/12, 7/12, 9/12, 
and 11/12 in Fe;W,. They are separated by 
distances of d=3.85A in Fe.2W and d=4.30A in 
Fe;W,. The distances of separation are de- 
termined by the positions of the intermediate W 
atoms. In both structures the remainder of the 
atoms just fill in the spaces between these layers, 
going into positions above and below x and y 
in Fig. 9, with most of the large tungsten atoms 
above and below position x where there is the 
more room. 

The fact that this network is common to the 
two structures indicates that it is structurally 
important and is therefore relatively stable. This 
conclusion is supported by the known fact that 
Fe;W, crystallizes in the form of plates parallel to 
the basal plane. Perhaps also significant is the 
Fe-Fe distance of 2.36A in the network as com- 
pared with the closest distance of approach of Fe 
atoms, 2.50A, in @ and y iron and other inter- 
metallic compounds. 

A corollary to the structural importance of 
these iron planes is that they will play an im- 
portant part in the formation of precipitate. 
These planes, from this point of view, would be 
the first thing to form from the matrix which 
would be characteristic of the precipitate. Once 
started, others would form parallel to the first 
and separated from it by distances determined 
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by the local tungsten concentration. Expected 
diffraction effects from this initial situation would 
be: ool maxima (using hexagonal indices analo- 
gous to the intermetallic compounds) which 
would be broad corresponding to variable inter- 
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Fic. 7. A plot of the positions and relative intensities of 
the diffraction patterns of Figs. 5 and 6, with the published 
data on Fe2W and Fe7W, for comparison. 


network spacing; and the more important hko 
lines, which would be sharp corresponding to 
rather definite atom spacings within the network. 
Since in this initial condition the Fe and W 
atoms between the networks will be randomly 
distributed among the possible x and y positions, 
hkl lines cannot be predicted. 

Figure 5 may be seen to correspond to this 
expected diffraction pattern. The sharp maxima 
at d=2.37A corresponds exactly to the common 
110 line of the intermetallic compounds. The 
broad band at d=3.94 to 4.40A may be con- 
sidered first order (001) diffraction from the 
networks spaced at varying distances corre- 
sponding to the spacings in FesW and Fe;W, of 
d= 3.85 and 4.30A, respectively. The broad band 
at d=2.00 to 2.18A is then second order (002) 
and is much more intense because other atoms, 
Fe and W, will lie just halfway between networks. 
This diffraction pattern, then, may be con- 
sidered to conform to what is expected of the 
initial stages of precipitation. 

In the course of time the initial condition will 
tend to assume the equilibrium condition which 
for the precipitate is crystalline Fe2W according 
to the equilibrium diagram. This process would 
require that the inter-network atoms diffuse 
around, the W atoms gradually occupying all the 
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x positions, leaving the y position for the Fe 
atoms. This process and the accompanying 
evening out of the composition would cause the 
inter-network spacing to assume a more definite 
uniform value tending toward 3.85A. The dif- 
fraction pattern should show this as a sharpening 
of the ool reflections and the appearance of hkl 
reflections because the W’s would be in definite 
position relative to the intra-network Fe’s. 

The diffraction pattern, Fig. 6, shows a stronger 
110 line at 2.37A, more definite 001 and 002 lines 
at 4.20 and 2.10A, respectively, and in addition, 
two broad lines which cannot be definitely 
assigned indices, but which correspond roughly to 
low order lines in either Fe2W or Fe;W,g (Fig. 7). 
Qualitatively, then, the above picture is con- 
firmed by the diffraction pattern of the inter- 
mediate stage. The 001 spacing, 4.20A seems 
larger than expected but the pattern of Fig. 7 is 
certainly not that of equilibrium material. 


V. CONCLUSIONS 


The lattice parameter measurements and the 
precipitate diffraction evidence indicate that age 
hardening in this system is caused by the pres- 
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Fic. 8. The structure of Fe2W and Fe7W, as determined by 
Westgren and Arnfelt (reference 4). 


ence in the hardened material of finely divided 
precipitate which is crystallographically distinct 
from the matrix. 

The precipitation reaction is of the discon- 
tinuous type, precipitation taking place in 
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restricted regions possibly associated with slip 
planes. 

A picture of the mechanism of precipitation, 
deduced from diffraction evidence and the struc- 
ture of the intermetallic compounds in the 
system, has been presented. 

Results of this work incidental to the main 
purpose are the lattice parameter—composition 
curve and the observation that the equilibrium 
solid solution composition at 700°C is probably 
1.4 atomic percent W. 

Since this picture of precipitation is different 
from those advanced for Al-Cu and Al-Ag, it may 
be worth while to summarize the possibly im- 
portant points of difference. The first con- 
sideration is that of atomic size. Cu is smaller 
than Al, Ag is almost exactly the same size, while 
W is considerably larger than Fe. In addition to 
affecting diffusion, the size factor plays an 
important part in the crystal structure of the 
intermetallic compounds and therefore in the 
nucleation mechanism. For instance, the segre- 
gation of large W atoms on any plane of the 
matrix containing a close packed direction would 
be unlikely. The second consideration is the 





structure of the matrix. Al-Cu and Al-Ag are 
face-centered cubic while Fe-W is body-centered 
cubic. The third possibly important factor is the 





COOOO0OCO 
Fic. 9. The distribution of Fe atoms in planes parallel 


to the basal plane of Fe2W and Fe;W¢. The unit cell in the 
plane at height z= { in Fe2W is outlined. 


electronic structure of the atoms. Fe and W are 
both transition elements and W has the highest 
binding energy of all the metals, while neither 
Al, Cu, nor Ag are transition metals. 

The author wishes to express his appreciation 
for the advice and criticism of Professor B. E. 
Warren under whose direction 
done. 
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